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PREFACE 


For  many  years,  I  have  been  fascinated  with  the 
profound,  yet  mysterious  capabilities  of  the  brain.  This 
marvelous  mass  of  tissue  is  able  to  perform  complicated  tasks 
involving  the  combination  of  sight,  touch,  hearing,  taste,  and 
smelling,  while  automatically  performing  life  sustaining 
functions.  Although,  these  tasks  are  often  complex,  many  of 
them  can  be  performed  by  young  children.  Yet,  how  the  brain  is 
able  to  perform  these  tasks  is  still  much  of  a  mystery. 

For  many  years,  researchers  in  the  area  of  brain  and 
central  nervous  system  research,  have  used  two  methods  to 
attempt  to  unlock  the  mysteries  of  the  brain.  The  first  method 
was  the  use  of  a  standard  electroenphalographic  ( EEG)  machine, 
using  scalp  electrodes,  to  measure  the  gross  electrical 
activity  of  the  brain  as  measured  on  the  surface  of  the  scalp. 
The  second  method  used  to  gain  some  insights  into  the 
operation  of  the  brain  was  the  meticulous,  and  painstaking 
measurements  of  the  responses  of  individuals  brain  cells  to 
specific  sensory  stimulus,  via  the  use  of  microelectrodes, 
a  extremely  slender  glass  or  metal  electrodes  which  are 
inserted  into  the  individual  cells. 

It  was  late  in  the  1960's  before  advancements  in 
cryogenics,  and  in  the  technology  of  superconducting  semicon¬ 
ductors  devices,  brought  about  a  new,  and  powerful  tool  to  be 
used  in  the  quest  for  understanding  of  the  brain.  This 
powerful  tool  is  the  Superconducting  Quantum  Interference 
Device  (SQUID)  Gradiometer.  This  device  allows  us  for  the 
first  time  to  measure  extremely  small  magnetic  field,  as  those 


associated  with  brain  activity.  I  am  grateful  to  have  the 
opportunity  to  work  with  such  a  state-of-the-art  research  tool 
in  my  research,  and  at  such  an  historical  time  as  soon  after 
the  successful  implantation  of  the  AFIT  array  integrated 
circuit,  the  first  cortical  implant  of  a  multiplexed 
multi-electrode. 

I  gratefully  acknowledge  the  support  of  Mr.  Gary  Reid 
of  the  Aerospace  Medical  Research  Laboratory  in  seeing  that 
equipment  was  available  to  me  as  needed.  Also,  my  thanks  go  to 
Mr.  Alva  Karl  and  the  Life  Science  personnel  of  the 
laboratory,  for  the  preparation  and  care  of  my  canine  subject, 
Ricky. 

My  sincere  appreciation  is  extended  to  many  people 
working  for  the  Systems  Research  Laboratories,  Inc.  Some  of 
there  employees  which  I  would  like  to  thank  especially  are  : 
Barbara  McFawn  for  her  overall  support  and  technical  advice, 
Karen  Peio  for  her  expert  instructions  in  use  of  laboratory 
equipment,  Donald  Stafford  for  his  advice,  and  maintenance  of 
the  SQUID  gradiometer,  Shelton  Unger  for  his  expert  software 
support.  Heather  Dippel  for  being  a  willing  MEG  subject, 
Regina  Palmer,  Melita  Ritter  for  their  support  in  my  data 
collection,  and  finally,  David  Holloway  for  his  graphics 
support. 

I  gratefully  recognize  Dean  L.  E.  Wolaver  for  his 
interest  and  work  as  a  committee  member,  and  his  instruction 
in  his  interesting  bioengineering  class. 

I  would  especially  like  to  thank  Dr.  Matthew  Kabrisky 
for  his  excellent  classes  on  the  central  nervous  system,  and 
in  sensory  communication.  These  classes,  along  with  my  other 


bioengineering  classes  motivated  me  to  quest  further  knowledge 
in  these  areas.  Also,  I  would  like  to  thank  him  for  support, 
and  encouragement  in  my  thesis  effort. 

Another  special  thanks  goes  to  my  thesis  advisor.  Dr. 
Charles  P.  Hatsell  for  his  expert  support,  guidance,  and 
technical  advice.  His  expert  knowledge  in  the  areas  of 
engineering  and  medicine,  and  his  high  degree  of  profession¬ 
alism,  sets  a  high  standard  for  a  researcher  to  follow. 

Finally,  I  would  like  to  express  my  sincere,  and 
equally  important  gratitude  to  my  wife.  Sherry  Murray,  and  our 
sons,  Robert  D.  Jr.,  and  John  N.  Murray,  for  their  patience 
and  understanding  during  the  long  work  hours,  both  at  home  and 
away  from  home. 


Robert  D.  Murray 
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ABSTRACT 


The  acquisition  of  a  Superconducting  Quantum 
Interference  Device  (SQUID)  gradiometer^  by— -the  Air  Force 
-AercrSp3jSe  Medical  Research  Laboratory,  along  with  the  previous 
research  oewducted  at  the  Air  Force  Institute  of  Technology 
which  led  to  the  first  cortical  implant  of  a  multiplexed 
multi-electrode  semiconductor  brain  electrode,  motivated  the 
use  of  the  SQUID  second-order  gradiometer  to  measure  visual 
evoked  responses  (VER)  from  humans,  and  from  the  canine  which 
had  previously  had  the  AFIT  brain  electrode  implanted  and 
removed.  Some  discussion  of  the  current  dipole  model  for  the 
fields  measured  from  the  brain  is  presented.  The  human  visual 
evoked  responses  picked  up  by  gradiometer  were  compared  to 
previous  human  electroencephalograph  visual  evoked  responses, 
and  were  found  to  have  similar  latencies.  The  canine  visual 
evoked  responses  which  were  measured  by  the  gradiometer  had 
similar  initial  waveform  and  latencies,  but  then  dampen  out 
more  rapidly  than  the  waveforms  measured  by  the  AFIT  brain 
chip.  Conclusions  about  the  use  of  a  SQUID  gradiometer  in 
measuring  visual  evoked  reponses,  along  with  comparisons  been 
the  EEG  and  AFIT  brain  chip  data  versus  the  SQUID  data  is 
discussed. KRecommendations  for  further  research  with  the  SQUID 


gradiometer  by  itself,  and  it  in  conjuction  with  the  AFIT 
brain  chip  are  discussed. 


CHAPTER  I 


INTRODUCTION 

BACKGROUND 

Ever  since  the  discovery  of  the  electrical  properties  of 
nerves  in  the  eighteenth  century,  scientists  have  been 
studying  this  phenomenon  in  the  hope  that  with  its  under¬ 
standing,  they  would  be  better  able  to  understand  how  the 
central  nervous  system  functions  (Ref.  1:2).  This  quest  to 
understand  the  central  nervous  system,  and  especial  y  the 
brain,  gave  birth  to  the  modern  science  of  neurology,  which  is 
the  branch  of  medicine  which  deals  with  the  nervous  system  and 
its  diseases,  and  to  the  modern  science  of  neurophysiology, 
which  is  the  branch  of  medicine  which  deals  with  the 
physiology  of  the  nervous  structure  of  the  body.  Even  though, 
through  the  years,  many  new  tools  were  developed  for  neural 
researchers,  these  tools  were  usually  improved  versions  of 
their  predecessors,  rather  than  new  innovations. 

Ever  since  Hans  Berger  discovered  human  brain  waves  in 
1924,  after  he  placed  small  metal  strips  on  the  scalps  of  his 
subjects  and  ran  wires  from  these  strips  to  a  sensitive 
galvanometer,  the  electroencephalograph  (EEG)  has  been  a  major 
tool  used  by  researchers  to  study  brain  activity  (Ref.  1:98). 
Many  improvements  have  been  made  to  the  standard  electro¬ 
encephalograph  over  the  years,  but  still  the  basic  concept  of 
the  EEG  is  measurement  of  electrical  potentials  on  the 
subjects  scalp  via  "gross  electrodes"  placed  at  specific 
locations  on  the  scalp,  and  then  interpreting  these 


measurements  in  the  hope  that  detailed  information  about 
mental  processes,  and  the  mental  and  physiological  condition 
of  the  subject  might  be  realized.  A  good  example  of  the  use  of 
the  EEG  as  a  diagnostic  tool  is  shown  in  Figure  1-1,  which 
shows  an  EEG  record  taken  during  a  brief  epileptic  fit  —  note 
the  spiky  waves  during  the  attack  (  starting  at  the  arrow 
marked  a  )  and  the  return  of  the  alpha  rhythm  after  eye 
closure  (  arrow  marked  b  )"  (Ref.  2:148).  For  more  background 
in  the  use  of  EEG,  see  references  3, 4, 5, 6, 7,  and  8. 
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FIGURE  1-1 

EEG  RECORD  OF  BRIEF  EPILEPTIC  FIT 
(  Ref.  2:148  ) 

Another  tool  which  researchers  have  used  to  study  brain 
activity,  is  the  microelectrode.  The  basic  microelectrode  is 
usually  a  extremely  slender  glass  or  metal  electrode  which  is 
inserted  into  or  near  an  individual  neuron.  Once  again,  it  is 
electrical  potentials  which  are  measured,  but  this  time  the 
firing  rates  of  action  potentials  of  the  neurons  are  analyzed. 
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This  tool  has  been  used  in  many,  and  varied  research  in  neuron 
activity.  Figure  1-2  shows  nerve  fiber  activity  signals  which 
were  taken  by  a  microelectrode  inserted  into  the  hand  of  a 
subject  lying  at  rest  (Ref.  9:141). 
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Nerve  fiber  activity  underlies  breathing,  blood  flow  and  skin  responses  to  disturbances  of 
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FIGURE  1-2 

EXAMPLE  OF  NERVE  FIBER  ACTIVITY 
TAKEN  WITH  A  MICROELECTRODE 
(  Ref.  9:141  ) 

For  more  background  in  the  applications  of  microelec¬ 
trodes  in  neural  research  see  references  9  thru  16. 

Although,  both  the  use  of  the  electroencephalograph  and 
the  use  of  the  microelectrode  has  provided  us  with  valuable 
information  in  the  study  of  the  central  nervous  system,  and 
related  electrical  phenomenon,  two  newer  innovations  have 
provided  additional  understanding  of  neural  activity.  These 
two  advance  measurement  devices  are  AFIT  Cortical  Implant  and 


the  Superconducting  Quantum  Interference  Device  (SQUID) 
Gradiometer . 


PIGURE  1-3 

AFIT  MULTI ELECTRODE  ARRAY 
(REF  18:55) 

In  October,  1982,  the  first  cortical  implant  of  a 
multiplexed  multi-electrode  semiconductor  brain  electrode  was 
successfully  implanted  on  the  brain  of  a  laboratory  beagle 
(Canis  familiaris)  (Ref.  17  and  Ref.  18).  It  is  an  x-Y 
addressable,  16  electrode  (4x4  array  )  multiplexed 


switching  electrode  based  on  contemporary  silicon  computer 
chip  technology  (Ref  17).  The  circuit  diagram  for  the 
electrode  is  shown  in  Figure  1-3. 

The  development  of  the  AFIT  implant  is  important,  since 
it  allows  us  to  measure  potentials  directly  from  a  very  small 
area  of  the  cortex  of  the  brain;  the  electrodes  were  spaced  30 
microns  apart  in  one  direction  and  100  microns  in  the 
orthogonal  direction. 

Now  that  we  have  a  device  which  can  measure  potentials 
directly  from  the  surface  of  the  cortex  of  the  brain,  we  can 
look  at  other  approaches  to  the  measurement  of  brain  activity. 
The  relatively  new  approach  to  the  measurement  of  brain 
activity  is  the  use  of  a  SQUID  gradiometer  to  measure  the 
magnetic  fields  emanating  from  the  brain.  This  approach  was 
slow  in  coming  because  we  didn't  have  the  technological 
advancements  necessary  to  take  these  measurements  until 
relatively  recent  years.  These  advancements  were  made  in  the 
fields  of  cryogenics,  and  the  very  low  temperature  supercon¬ 
ducting  semiconductor  physics,  which  developed  Josephson 
Junction  electronic  devices  to  be  used  in  the  Superconducting 
Quantum  Interference  Device.  A  discussion  of  the  theory 
involved  in  the  operation  of  a  SQUID  gradiometer  will  be  later 
(Chapter  2). 

Having  now,  both  the  AFIT  cortical  implant  electrode 
research  continuing,  and  research  at  the  Air  Force  Aerospace 
Medical  Research  Laboratory  with  a  second-order  SQUID 
gradiometer,  I  decided  to  take  a  series  of  magnetic  measure¬ 
ments  with  the  SQUID  gradiometer,  and  compare  the  results  with 


predicted  values,  and  with  measurements  taken  by  the  AFIT 
cortical  implant  electrode. 

STATEMENT  OF  THE  PROBLEM 

The  primary  emphasis  of  this  research  is  the  comparison 
of  measurements  taking  with  a  second-order  SQUID  gradiometer 
to  waveforms  either  predicted  mathematically,  or  to  waveforms 
measured  in  the  past  with  an  EEG  or  a  AFIT  cortical  implant 
electrode. 

The  first  area  of  investigation  was  comparing  the 
predicted  mathematical  magnetic  and  potential  fields  waveform 
emanating  from  a  current  dipole  submerged  in  a  saline 
solution. 

Next,  magnetoencephalographic  visual  evoked  responses 
were  measured  from  several  human  subjects,  and  compared  with 
the  basic  EEG  visual  evoked  response  waveform. 

Finally,  magnetoencephalographic  visual  evoked 
responses  were  measured  from  the  same  laboratory  beagle  which 
had  had  electroencephalographic  visual  evoked  responses 
measured  by  the  AFIT  cortical  implant  electrode. 

SCOPE 

When  the  dewar  which  held  the  liquid  helium  for  the 
SQUID  gradiometer  failed  to  maintain  superconducting  tempera¬ 
tures,  the  ensuing  downtime,  along  with  the  growing  interest 
to  measure  magnetoencephalographic  visual  evoked  responses 
from  the  same  laboratory  beagle  which  had  had  electroenceph¬ 
alographic  visual  evoked  responses  measured  by  the  AFIT 


cortical  implant  electrode,  I  decided  to  place  my  primary 
emphasis  in  taking  measurements  with  the  SQUID  gradiometer.  I 
especially  placed  my  emphasis  on  measuring  visual  evoked 
responses  from  the  laboratory  beagle,  instead  of  doing  a 
feasibility  study  of  simultaneous  or  near  simultaneous 
electroencephalographic  and  magnetoencephalographic  measure¬ 
ments  from  a  brain. 

The  problem  with  the  SQUID  gradiometer  downtime  also 
reduced  the  amount  of  time  I  had  to  collect  data,  but  suf¬ 
ficient  data  was  collected  to  make  conclusions  possible. 

APPROACH  AND  PRESENTATION 

This  research  is  presented  in  essentially  chronological 
order,  with  each  chapter  representing  a  phase  of  research. 

Chapter  II  gives  a  brief  discussion  into  the  theory 
behind  the  SQUID  gradiometer,  along  with  some  of  it  proposed 
applications. 

Chapter  III  gives  a  look  at  the  basic  current  dipole 
model  for  the  brain,  along  with  the  predicted  and  measured 
results  for  a  current  dipole  submerged  in  a  saline  solution. 

Next,  Chapter  IV  discusses  some  basic  theory  behind 
visual  evoked  response  measurements,  along  with  human  magneto¬ 
encephalographic  visual  evoked  responses  waveforms  compared  to 
the  standard  BEG  waveform  for  a  visual  evoked  response. 

Chapter  V  details  the  magnetoencephalographic  measure¬ 
ment  of  the  laboratory  beagle's  visual  evoked  response,  along 
with  the  comparison  of  the  magnetoencephalographic  visual 
evoked  response  to  the  electroencephalographic  visual  evoked 


response. 

Chapter  VI  draws  conclusions  about  the  use  of  the  SQUID 
gradiometer  to  measure  magnetoencephalographic  visual  evoked 
responses,  and  summarizes  the  results  of  this  research. 

Finally,  Chapter  VII  makes  recommendations  for  further 
experiments  using  the  SQUID  gradiometer. 


CHAPTER  II 


SQUID  G RADIOMETER  THEORY 


THEORY 

SQUID  is  an  acronym  of  "superconducting  quantum 
interference  device."  In  this  chapter  I  will  first  discuss 
general  SQUID  theory,  and  then  I  will  discuss  some  parameters 
and  details  of  the  specific  SQUID  gradiometer  which  I  used  in 
measuring  the  visual  evoked  responses. 

First  lets  look  at  the  SQUID  gradiometer  as  a  black  box 

to  see  how  it  behaves  when  magnetic  flux  crosses  perpendicular 

to  the  input  detection  coil.  The  input  to  this  black  box  is 

the  current  which  the  external  magnetic  flux  cutting  the  input 

detection  coil  induced  into  the  detection  loop.  The  output  of 

the  black  box  is  a  voltage  which  is  proportional  to  the  amount 

of  magnetic  flux  which  has  perpendicularly  penetrated  the 

effective  area  of  the  input  loops.  So,  we  can  say  that  it 

simply  behaves  like  a  cur rent-to-voltage  converter.  It  offers 

very  high  gain,  extremely  low  noise,  wide  dynamic  range, 

excellent  linearity,  and  a  wide  bandwidth,  extending  down  to 

dc.  The  current-to-voltage  conversion  factor  for  a  SQUID 

7 

gradiometer  is  usually  approximately  10  V/A  (Ref.  26:70). 

Well,  we  see  from  the  black  box  outlook,  the  SQUID 
gradiometer  can  give  us  a  voltage  which  is  proportional  to  a 
small  magnetic  flux  detected  by  its  input  detection  coil.  This 
is  an  interesting  conclusion,  but  what  is  the  significance  of 
the  SQUID  gradiometer  which  makes  it  the  major  tool  in  the 
area  of  biomagnetic  measurement?  It  is  the  SQUID  gradiometer ' s 
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capability  to  detect  a  magnetic  field  as  weak  as  10 
which  makes  it  so  attractive  in  many  areas  of  research.  That 

weak  magnetic  field  is  over  10^  times  smaller  than  the 

magnetic  field  of  the  earth. 

Now  that  we  have  an  idea  of  some  of  the  capabilities  of 
a  SQUID  gradiometer,  lets  discuss  in  greater  detail  the 

attributes  of  a  SQUID  gradiometer  which  gives  it  those 
capabilities.  First,  let  us  start  by  disecting  the  acronym 
"SQUID*  into  its  components.  The  "S"  stands  for 

superconducting.  This  device  uses  a  special  semiconductor 
device  which  operates  at  superconducting  temperatures,  that 
is,  temperatures  usually  below  20°  Kelvin.  This  special  device 
is  called  a  superconducting  insulator-barrier  junction,  or 
more  commonly  as  a  Josephson  junction,  after  the 
Noble-prize-winning  English  gratuate  student,  Brian  Josephson, 
who  predicted  the  junctions  properties  in  1962  (Ref  27:42).  In 
the  literature,  the  Josephson  junction  is  referred  to  as  a 
"weak  link"  device.  Josephson  junctions  which  we  find  in  SQUID 
gradiometers  many  times  rely  on  the  oxide  at  the  surface  of  a 
sharp  niobium  point  where  it  presses  against  a  flat  niobium 
surface.  They  are  called  "point  contact"  Josephson  junctions. 
The  superconducting  transition  temperature  for  niobium,  that 
is,  the  temperature  at  which  niobium  resistance  drops 
dramatically  to  approach  zero  resistance,  is  9.2°  Kelvin.  The 
resistance  of  the  superconducting  material  is  not  exactly 
zero,  since  you  still  have  some  molecular  motion  even  at  9.2° 
Kelvin.  To  get  zero  resistance,  the  superconducting  material 
must  be  at  0°  Kelvin.  Even  though  the  resistance  is  not 


exactly  equal  to  zero  for  superconductors  operating  at  9.2° 

Kelvin,  by  using  rings  with  very  low  inductance  and  a  very 

long  measurement  time,  together  with  a  high  sensitivity 

current  detector,  it  has  been  possible  to  demonstrate  that  the 
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residual  resistivity  of  a  superconductor  is  at  least  10 
smaller  than  the  resistivity  of  copper  at  room  temperaure. 
Therefore,  for  all  practical  purposes,  we  can  consider  that 
the  superconducting  material  has  zero  resistances  (Ref. 
26:72).  To  insure  that  the  SQUID  remains  superconducting,  it 
is  usually  immerged  in  a  bath  of  liquid  helium.  The  4.2 
degrees  Kelvin  temperature  of  the  liquid  helium  bath  is 
adequate  to  keep  the  SQUID  superconducting  (Ref.  26:69). 

But  what  effects  of  this  superconducting  Josephson 
junction  allows  it  to  be  used  in  a  SQUID  gradiometer?  First 
lets  review  the  Meissner  effect  of  a  superconductor.  Back  in 
1933  in  Berlin,  two  men,  Meissner  and  Ochsenfeld,  found  that 
the  interior  of  a  superconducting  material  is  free  of  magnetic 
flux  independent  of  whether  the  magneic  field  is  applied  above 
or  below  Tc,  the  superconducting  transition  temperature.  This 
property  is  now  commonly  known  as  the  Manner  effect.  In  the 
presence  of  a  magnetic  field,  a  current  in  a  very  thin  surface 
layer  is  generated  in  such  a  way  that  its  field  compensates 
the  field  inside.  This  is  distinctly  different  that  for  an 
ideal  case  of  Lenz's  law,  since  the  expulsion  of  flux  from  the 
Meissner  effect  occurs  even  in  a  constant  field  as  a  sample  is 
cooled  through  its  transition  temperature.  The  superconductor 
behaves  as  an  ideal  diamagnetic  material  with  susceptibility 
equal  to  -1  (Ref.  26:72). 


The  "QU"  and  the  "I"  in  the  acronym  comes  once  again 
from  the  effects  of  a  Josephson  junction  operating  at 
superconducting  temperatures.  The  "QU",  as  mention  above, 
stands  for  quantum,  and  the  "I"  stands  for  interference.  And 
finally,  the  "D"  just  stand  for  device.  Thus  the  name 
■superconducting  quantum  interference  device"  derives  from  the 
fact  that  quantum  interference  phenomena  can  be  observed  with 
SQUIDs  on  a  macroscopic  scale.  For  example,  the  current 
through  a  two-junction  dc-SQUID  varies  periodically  in  a 
manner  which  is  analogous  to  the  optical  interference  pattern 
from  two  slits.  At  the  Josephson  junction  insulating  barrier, 
the  amplitude  of  the  "wave  function"  describing  the 
wavelike-properties  of  the  superconducting  electrons  cannot 
drop  abruptly  to  zero,  because  then  according  to  quantum 
theory  the  electrons  would  have  an  enormously  high  kinetic 
energy,  which  is  not  permitted.  Instead,  the  wave  function 
penetrates  into  the  wall  a  short  distance  as  it  decays 
smoothly  toward  zero.  However,  if  the  barrier  is  sufficiently 
thin,  with  superconducting  electrons  on  both  sides,  the 
decaying  wave  functions  from  the  two  sides  overlap 
appreciable,  and  this  overlap  implies  that  superconducting 
electrons  can  flow  through  the  insulator  without  any 
resistance.  This  effect  is  known  as  tunnelling  because  the 
electrons  penetrate  through  a  potential  energy  barrier,  where 
by  the  classical  laws  of  physics  they  are  forbidden  to  go 
(Ref.  26:74).  But  how  did  the  electrons  manage  to  penetrate 
through  the  potential  energy  barrier?  Well,  this  phenomena  can 
be  explained  by  the  microscopic  theory  of  the  superconducting 


state  which  was  published  by  Cooper  and  Shrieffr  in  1957.  The 
basis  for  this  theory  is  the  consideration  that  at  very  low 
temperature  a  weak  attractive  force  between  electrons  becomes 
more  important  than  the  disordering  effect  of  thermal 
agitation,  and  they  bind  together  in  pairs  called  Cooper 
pairs.  It  is  these  Cooper  pairs  which  can  tunnel  through  the 
barrier  as  a  unit.  The  Cooper  pairs  are  not  scattered  in  a 
superconductor  whereby  wave  coherence  is  conserved  over 
infinitely  long  distances.  Interference  effects  can  thus  be 
observed  on  a  macroscopic  scale.  It  can  be  shown  that  the 
magnetic  flux  (D^  passing  through  the  area  of  the  weak  link  can 
be  taken  in  account  when  calculating  Jc,  the  critical  current 
density  for  the  junction.  Jc  becomes  the  periodic  function 
shown  in  Equation  2-1. 


Jc  "  Jc(0)  * 


I  sin(Pl  *  <Dj/0o)  I 


PI  *  (Dj/<t>0 


where 

Jc  is  the  critical  current  density 
PI  is  3.1415926 

is  the  magnetic  flux  passing 
through  the  weak  link 
<D0  is  the  flux  quantum  which  is  equal 
to  2.07  *  lO'15  Weber 


This  formula  has  the  same  form  as  the  optical 
Fraunhofer  diffraction  from  a  slit.  Because  (JK  is  much  smaller 
than  <D,  the  magnetic  flux  threading  through  the  superconducting 
loop,  the  interference  pattern  from  a  SQUID  is  modulated  by 
"diffraction"  occurring  at  the  weak  link. 

If  wires  are  connected  to  each  side  of  the  weak  link, 
the  voltage  measured  will  be  proportional  to  the  externally 
applied  magnetic  flux. 
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FIGURE  2-1 

DC  SQUID  MAGNETOMETER 

Now  that  we  have  discussed  how  a  Josephson  junction  can 
be  used  as  a  weak  link  in  a  superconducting  ring  so  as  to 


detect  small  magnetic  fields,  lets  talk  about  the  rest  of  a 
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SQUID  gradiometer.  There  are  two  types  of  SQUID  gradiometers . 
The  first  type  of  SQUID  gradiometer  which  I  will  discuss  is 
the  DC  SQUID  gradiometer.  Figure  2-1  shows  a  diagram  of  a  DC 
SQUID  gradiometer. 

The  DC  SQUID  gradiometer  consist  of  a  superconducting 
ring  interrupted  by  two  Josephson  junctions,  the  name  follows 
from  the  fact  that  in  operating  this  two-junction  SQUID  only  a 
dc  bias  current  is  necessary.  The  voltage  is  taken  off  of  the 


superconducting  ring  and  passed  through  a  step-up  transformer 
to  the  lock  in  amplifier  and  then  amplified  by  a  dc  amplifier. 
Part  of  the  output  is  sent  by  to  a  feedback  coil  too  provide 
the  feedback  necessary  to  allow  the  lock-in  amplifier  to 
maintain  flux-locked  mode.  Flux-lock  mode  and  the  lock-in 
amplifier  will  be  covered  later. 
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FIGURE  2-2 


RF  SQUID  MAGNETOMETER 


The  second  type  of  SQUID  gradiometer  is  the  RF  SQUID 
gradiometer.  A  diagram  of  a  RF  SQUID  gradiometer  is  shown  in 
Figure  2-2.  Once  again  a  SQUID  loop  is  used,  but  this  time 
only  one  Josephson  junction  is  used.  The  RF  SQUID  gradiometer 
is  the  type  usually  commercially  produced.  In  this  case,  the 
SQUID  loop  is  weakly  coupled  to  a  resonant  tank  circuit  tuned 
to  a  sinusoidal  pump  current  in  the  radio  frequency  range.  The 
rf  current  source  supplies  a  "pumping*  current  to  a  resonant 
tank  circuit.  The  current  flowing  between  the  inductor  and 
capacitor  of  the  tank  circuit  exceeds  the  pumping  current  by 
the  quality  factor  "Q"  of  the  tank  circuit.  This  circulating 
current  generates  an  rf  voltage  across  the  parallel 
combination  of  the  inductor  and  capacitor  which  is  usually 
enhanced  by  a  preamplifier,  rectified  by  an  rf  detector,  and 
then  smoothed  by  a  low-pass  filter  (usually  a  resistor  and 
capacitor)  to  produce  an  average  voltage  that  is  proportional 
to  the  amplitude  of  the  tank  circuit's  rf  voltage.  Note  that 
the  voltage  is  affected  by  the  state  of  the  SQUID,  because  the 
inductor  of  the  tank  circuit  is  magnetically  coupled  to  the 
SQUID.  The  changing  field  of  the  inductor  induces  a  current 
that  flows  around  hte  SQUID  loop,  and  this  rf  current  is 
superimposed  on  a  second  current  representing  the  input  signal 
that  results  from  the  field  applied  by  the  input  coil.  In  the 
simplest  case,  the  non-linear  response  of  the  junction  has  the 
effect  of  changing  the  inductance  of  the  SQUID  due  to  the 
presence  of  this  signal  current.  This  change  in  SQUID 
inductance  slightly  alters  the  inductance  in  the  tank  circuit 
to  which  it  is  magnetically  coupled,  and  thereby  it  affects 


the  tank  voltage.  This  xs  a  parametric  effect  (Ref.  26:78). 

Note  that  the  diagrams  in  Figure  2-1  and  Figure  2-2  are 
magnetometers,  and  not  gradiometers .  For  a  magnetometer  to 
become  a  gradiometer,  it  must  have  a  specific  coil 
configuration.  Figure  2-3  shows  four  coil  configurations. 


o)  b)  c) 


FIGURE  2-3 
COIL  CONFIGURATIONS 

The  coil  in  Figure  2-3(a)  is  a  magnetometer  coil  for 
measuring  the  average  flux  density  over  the  area  of  the  coil. 
Both  coils  ih  Figure  2—3 ( b )  are  used  in  first  order 
gradiometers,  that  is  ,  with  gradiometers  which  will  measure 
flux  differences  of  approximately  one  component  of  the  field 
gradient.  They  allow  measurements  of  local  biomagnetic  fields 
in  the  presence  of  much  stronger  but  uniform  fields  from 
distant  sources,  since  for  distant  sources  the  change  in  the 


partial  derivative  of  the  magnetic  field  with  respect  to 
distance  in  the  "z"  direction  is  approximately  equal  to  zero. 
The  last  coil  configuration  which  is  shown  in  Figure  2-3  is 
the  second-order  gradiometer  configuration  which  is  better 
suited  to  measure  extremely  small  local  magnetic  fields  while 
in  the  presence  of  higher  noise  levels,  since  not  only  are  the 
uniform  fields  rejected,  but  also  the  first-order  fields  are 
rejected. 

Because  of  its  ability  to  select  specific  magnetic 
fields  out  of  a  noisy  background,  the  second-order  SQUID 
gradiometer  is  normally  used  to  taxe  biomagnetic  measurements, 
including  all  the  measurements  reported  in  this  report. 

The  second-order  SQUID  gradiometer  can  also  be  used  to 
analyze  heart  activity  by  producing  a  magnetocardiogram  (MCG) , 
analyze  muscle  activity  by  producing  a  magnetomyogram  ( MMG ) , 
and  monitoring  movement  of  the  eye  by  producing  a 
magnetooculogram  (MOG). 

Before  leaving  this  chapter,  I  will  discuss  the 
operation  of  the  SQUID  gradiometer  in  "Flux-lock  mode",  and 
some  information  on  the  particular  SQUID  gradiometer  which  I 
used  to  take  measurements  for  this  report.  To  use  the  SQUID  in 
practical  instruments,  a  linearization  circuit  is  necessary  to 
minimize  the  nonlinearity  and  also  maximize  the  dynamic  range 
by  avoiding  the  periodicity  of  the  SQUID  responce.  For  this 
reason,  the  SQUID  normally  is  used  as  a  null  detector  in  a 
feedback  loop.  Usually,  an  audio  frequency  generator 
introduces  a  modulation  flux  into  the  SQUID.  The  amplitude  of 
the  modulation  is  approximately  one-fourth  of  a  flux  quantum. 


The  choice  of  the  modulation  frequency  depends  on  the  desired 
bandwidth  for  a  measurement,  because  this  frequency  has  to  be 
higher  than  the  highest  frequency  component  of  the  signal  to 
be  measured.  Usually  it  is  chosen  to  be  about  50  kHz  or 
higher.  A  phase-sensitive  detector  (PSD)  extracts  the  portion 
of  the  SQUID  electronics  output  at  this  modulation  frequency 
and  rectifies  it;  and  this  rectified  output  of  the  PSD  goes  to 
an  integrator  (low-pass  filter).  The  PSD  and  integrator 
comprise  what  is  known  as  a  "lock-in"  amplifier.  Figure  2-4 
showns  a  block  diagram  of  a  lock-in  amplifier  (Ref.  28:597). 


FIGURE  2-4 


DIAGRAM  OF  A  LOCK-IN  AMPLIFIER 


Finally,  I  will  discuss  the  S.H.E.  Corporation  System 


330X  Second-Order  SQUID  gradiometer,  which  I  used  in 
collecting  data  for  this  report. 

As  shown  in  Figure  2-5,  the  Model  330  RF  Unit  uses  a 
single  weak  link  made  from  niobium.  Its  rf  coil  requires  one 
microamp  to  change  flux  by  one  flux  quantum.  Its  input  coil 
has  a  current  sensitivity  of  .1  microamp  per  flux  quantum.  It 
uses  a  RF  bias  frequency  of  19  MHz  with  a  50  kHz  square  wave 
modulation.  The  input  has  an  overall  gain  of  83  db. 


DIAGRAM  OF  S.H.E.  MODEL  300  RF  UNIT 


Its  bandwidth  is  plus  and  minus  150  kHz  at  carrier 
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CHAPTER  III 


DIPOLE  MODEL  WITH  TESTING 


BACKGROUND 


In  this  chapter,  I  will  first  discuss  the  dipole  theory 
as  the  source  of  electrical  activity  in  the  brain.  Next,  I 
will  discuss  the  single  dipole  model.  In  this  case,  it  has 
been  placed  in  a  bowl  immerged  in  a  saline  solution.  The 
predicted  values  for  the  magnetic  field  emanating  from  the 
dipole  are  plotted,  and  compared  with  the  plot  of  the  measure¬ 
ments  taken  by  the  second-order  SQUID  gradiometer.  Finally, 
the  surface  potential  predicted  values  are  plotted,  and 
compared  with  the  plot  of  the  potentials  measured  off  the 
surface  of  the  saline  solution,  by  using  an  oscilloscope  with 
input  differential  amplifiers. 


DIPOLE  MODEL  FOR  THE  BRAIN 

With  approximately  60  to  75  percent  of  cortical  neurons 
oriented  perpendicular  to  the  cortical  surface,  and  a  large 
number  of  interconnections,  it  then  appears  that  the  scalp 
potentials  are  due  to  current  dipole  layers  occupying  various 
surface  areas  of  the  cortex  (Ref.  5:273).  For  more  background 
in  reference  to  the  dipole  model  for  the  brain  see  references 
5,6,  and  21  thru  24. 


SINGLE  DIPOLE 


For  this  part  of  my  research,  I  immerged  a  .5  cm. 


current  dipole  in  a  solution  of  .464  percent  saline.  Figure 


3-1  shows  a  diagram  of  the  current  dipole  experiment  appa¬ 
ratus,  while  Figure  3-2  shows  a  picture  of  the  apparatus.  This 
apparatus  included  two  48  cm.  by  31  cm.  by  1  cm.  pieces  of 


Plexiglas,  an  approximately  3  liter  glass  bowl,  a  2  cm.  by  2 
cm.  by  3  cm.  block  of  lucite,  and  a  3  meter  piece  of  shielded 
wire.  I  made  a  1/2  cm.  current  dipole  on  one  end  of  the  wire, 
and  connected  the  other  end  of  the  wire  thru  a  1000  ohm 
resister  to  a  BNC  connector.  On  one  of  the  large  pieces  of 
plexiglass  I  etched  a  40  cm.  by  25  cm.  grid.  This  grid  was 
graduated  in  1  cm.  increments.  From  the  top  left-hand  corner 
of  this  grid,  I  marked  the  origin  for  the  cartesian  coordinate 
system  to  be  down  12  cm.  and  right  20  cm.  From  the  origin  of 
this  cartesian  coordinate  system,  I  chose  the  direction  to  the 
right  to  be  the  positive  "y*  direction,  and  the  direction  down 
to  be  the  positive  "x*  direction.  The  positive  "z"  direction 
was  in  the  direction  above  the  apparatus.  This  piece  of  plexi¬ 
glass  became  the  bottom  of  my  apparatus,  and  was  used  to 
correctly  align  the  apparatus  to  the  desired  measurement 
point.  For  example,  to  align  the  apparatus  to  the  "x,y"  coor¬ 
dinate  (3,-5),  the  point  on  the  grid  corresponding  to  the 
intersection  of  the  "x"  equal  to  positive  3  line  and  the  "y* 
equal  to  negative  5  line  is  placed  on  top  of  a  center  of 
measurement  point  which  was  initially  determined  before  the 
dipole  apparatus  was  placed  under  the  SQUID  gradiometer.  This 
initial  center  of  measurement  point  was  determined  by  placing 
the  SQUID  gradiometer  in  its  vertical  position  and  placing  a 
mark  at  the  point  on  the  measurement  table  at  which  the  line 
projecting  from  the  center  of  the  SQUID  gradiometer  detection 


The  bowl  was  filled  with  the  .464  percent  saline 
solution  and  placed  between  the  pieces  of  plexiglass.  Four 
brass  screws  with  wingnuts  were  used  to  fasten  the  bowl  in 
place.  The  center  of  this  glass  bowl  was  determined  and 
marked.  This  mark  was  placed  on  top  of  the  origin  mark  on  the 
bottom  piece  of  plexiglas  so  that  the  dipole  would  be  centered 
in  the  cartesian  coordinate  system. 

Before  the  apparatus  was  completely  built,  I  calculated 

mathematically  the  predicted  values  of  the  magnetic  field 

strength  (H)  using  the  equation  derived  by  Cuffin  and  Cohen  in 

there  paper,  "Magnetic  Fields  of  a  Dipole  in  Special  Volume 

Conductor  Shapes.  (Ref.  24:374)"  Equation  3-1  shows  the 

formula  for  H  . 

z  • 


X*Py 


4*PI* ( X*X 

where 

Py  =  dipole  moment  = 
PI  =  3.14159 
X  =  distance  in  the 
Y  =  distance  in  the 
Z  =  distance  in  the 
D  =  depth  dipole  is 


+  Y*Y  +(D+Z)*(D+Z) )3/2 

current  *  dipole  length 

X  direction 
Y  direction 
Z  direction 
below  surface 


For  my  experiment. 


Py  =  ( ,008)*( .005) 


.00004  and  the 


value  of  D  was  .015  meters.  To  calulate  the  magnetic  flux 
density  (B),  I  used  Equation  3-2. 


Bz  =  Po*H 

where 

Bz  =  magnetic  flux  density 
=  4*PI*10“7 

Hz  =  magnetic  field  strength 

EQUATION  3-2 


The  calculated  magnetic  flux  density  (B)  values  for  the 
plus  and  minus  10  cm.  grid  is  located  in  Appendix  A,  along 
with  Table  A-l  which  shows  a  subset  of  the  measurements  which 
were  taken,  and  their  associated  predicted  value.  As  I  was 
analyzing  the  calculated  values  versus  the  measured  data  I  was 
looking  for  an  isofield  contour  shape  to  correlate  to  each 
other.  Figure  3-3  shows  a  current  dipole  with  its  associated 
magnetic  field  B  and  volume  current  density  Jv  Figure  3-4 
shows  the  plot  of  the  calculated  magnetic  field  for  the 
dipole.  Note  that  the  calculated  magnetic  fields  isofield 
contours  forms  two  lobes.  One  lobe  is  pointed  in  the  positive 
"x"  direction,  and  the  other  lobe  is  pointed  in  the  negative 
"x"  direction. 
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FIGURE  3-3 
CURRENT  DIPOLE 

Also  note  that  the  positive  "y"  direction,  that  is,  the 
direction  which  the  dipole  is  pointing,  is  to  the  right,  and 
the  positive  "x"  direction  is  toward  the  bottom  of  the  plot. 
This  plot  was  then  compared  to  the  plot  obtained  by  plotting 
the  data  measured  by  the  SQUID  gradiometer.  This  plot  is  shown 
in  Figure  3-5.  Although,  the  plot  differs  in  absolute  isofield 
value,  and  the  wave  shape  is  distorted,  we  can  still  see  that 
it  is  still  directly  correlated  to  the  calculated  field  plot. 
There  are  several  factors  involved  in  the  distortion  of  the 
isofield  contour  lines  in  the  measured  magnetic  field  plot. 
First,  we  must  remember  that  the  calculated  isofield  contour 
plot  did  not  take  into  account  the  edge  effects  which  are 
introduced  into  the  measurements  since  the  bowl  sides 
restricts  the  medium  to  a  finite  distance  in  the  "x"  and  "y" 
directions,  even  though  the  calculations  are  based  on  an 
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Also,  the  lower  lobe  is  greatly  distorted,  as  expected, 
because  of  the  lucite  block  which  is  used  to  mount  the  dipole 
is  at  the  negative  side  of  the  dipole,  so  the  volume  are 
forced  to  go  around  this  block  and  in  doing  so  get  distorted. 
The  difference  in  the  absolute  value  of  the  calculated  mag¬ 
netic  field  versus  the  measured  magnetic  field  is  partially 
due  to  the  edge  effects  introduced  into  the  measurements,  and 
also  the  possible  inaccuracies  in  the  positioning  of  the 
apparatus  with  respect  to  the  SQUID  gradiometer. 

Since  later  in  this  report  I  will  be  comparing  magneto- 
encephalographic  measurements  of  visual  evoked  responses  to 
electroencephalographic  measurements  of  visual  evoked 
responses,  my  single  dipole  experiment  would  not  be  complete 
without  predicting  the  surface  electrical  potential  of  the 
saline  solution,  and  then  compare  the  predicted  values  to  the 
values  I  obtain  from  measuring  the  surface  potential  using  a 
Tektronix  7A22  Differential  Amplifier  plugged  into  a  Tektronix 
7623A  Oscilloscope.  As  for  all  of  the  experiments  with  the 
single  dipole  apparatus,  the  Textronix  503  Function  Generater 
supplied  the  8  milliamps  of  current  to  the  dipole. 

In  order  to  solve  for  the  surface  potential  caused  by  a 
dipole  which  is  submerged  a  depth  "d"  below  the  surface  of  the 
solution,  let  us  first  consider  the  case  shown  in  Figure  3-6. 
In  Figure  3-6,  you  see  a  single  dipole  of  length  2  "e*.  The 
dipole  is  centered  at  the  origin,  but  pointing  in  the  "y" 
direction.  For  this  case,  the  conductivity  is  the  same 
throughout  the  medium.  Once  we  choose  an  arbitrary  point,  say 
P,  we  can  apply  the  potential  function  for  a  dipole  to 


determine  the  potential  at  point  p  due  to  the  dipole. 
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FIGURE  3-6 
SINGLE  DIPOLE  MODEL 


The  potential  at  point  P  due  to  the  dipole  can  be 
written  in  the  form  shown  in  Equation  3-3: 


p/e 

<t>d  =  -  * 

4*PI*SIGMA 


1 


EQUATION  3-3 


where 


p  =  dipole  moment 
e  =  half  the  length  of  the  dipole 


PI  =  3.14159 


R+  =  (  X*X+(Y-e)*(Y-e)  +  Z*Z)1//2 
R_  =  {  X*X+(Y+e)*(Y+e)  +  Z*Z)1//2 

=  potential  at  point  P  due  to  the  dipole 
SIGMA  =  the  conductivity  of  the  medium 

If  we  take  the  limit  of  the  expression  above  as  "e" 
goes  to  zero,  and  do  some  algebra,  we  get  the  expression  shown 
below  in  Equation  3-4: 


2*P*Y 


4*PI*SIGMA*R3 

EQUATION  3-4 

So  now  we  have  an  equation  for  the  potential  at  an 
arbitrary  point  P  due  to  a  single  dipole  in  a  infinite  medium. 
Now  we  can  use  the  above  equation,  and  the  superposition 
principle  to  develop  the  equation  for  the  potential  at  an 
arbitrary  point  P  due  to  two  dipoles.  Figure  3-7  shows  the  new 
problem  we  are  working,  that  is,  two  dipoles  are  separated  by 
2  "d",  but  they  are  orientated  in  the  same  direction.  Using 
the  superposition  principle,  we  get  Equation  3-5,  which  is 
shown  below,  as  the  solution  to  the  problem  of  two  dipoles 
placed  2  "d"  apart  and  with  the  same  orientation. 

The  calculated  potentials  are  located  in  Appendix  B, 
along  with  Table  B-l,  which  contains  a  subset  of  the  measure¬ 
ments  which  were  taken,  and  their  associated  predicted  value. 
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EQUATION  3-5 


where 

p  =  dipole  moment 
PI  =  3.14159 

R+  =  (  X*X+Y*Y+(Z-D)*<Z-D) )3/2 

R_  ■  (  X*X+Y*Y+{Z+D)*(Z+D) )3/2 

<t>d  =  potential  at  point  P  due  to  the  dipoles 

SIGMA  =  conductivity  of  the  medium 

Once  again,  we  will  analyze  the  plot  of  the  calculated 
dipole  surface  potential  shown  in  Figure  3-8,  along  with  the 
plot  of  the  measured  dipole  surface  potential  shown  in  Figure 
3-9.  Once  again,  the  positive  "y"  direction  is  toward  the 
right  of  the  plot,  and  the  positive  "x"  direction  is  toward 
the  bottom  of  the  plot.  In  Figure  3-8,  we  see  two  lobes  as  we 
saw  in  the  calculated  magnetic  field  plot  of  Figure  3-4,  but 
this  time  the  lobes  are  pointing  in  the  plus  and  minus  "y" 
direction  instead  of  the  plus  and  minus  "x"  direction.  Once 
again,  we  see  that  the  lobe  associated  with  the  positive  "y" 
direction  is  present,  but  that  the  lucite  block  effects  is 
once  again  present. 


CALCULATED  DIPOLE  SURFACE  POTENTIAL 


CALCULATED  DIPOLE  SURFACE  POTENTIAL 


Some  of  the  inaccuracies  in  both  the  wave  shape  and 
the  absolute  value  of  the  measured  data  versus  the  calculated 
values  is  due  to,  once  again,  the  edge  effects  of  having  a 
finite  volume,  rather  than  a  semi-infinite  volume.  Also, 
positioning  of  the  measurement  probe  may  have  had  some 
affects.  And  finally,  as  shown  in  Equation  3-5,  some  inac¬ 
curacies  may  have  been  introduced  by  using  an  estimated  value 
for  SIGMA,  the  conductivity  of  the  saline  solution. 


CHAPTER  IV 


HUMAN  VISUAL  EVOKED  RESPONSE 

BACKGROUND 

Human  brain  waves  were  discovered  in  1924  by  Hans 
Berger  of  the  University  of  Jena,  Ausria  (Ref.  1:98).  Since 
then  many  techniques  have  been  developed  to  study  brain 
activity.  One  major  advancement  in  the  area  of  brain  research 
was  the  development  of  the  electroencephalograph.  Through  the 
years,  the  electroencephalograph  has  had  many  applications, 
including  its  use  to  measure  visual  evoked  responses.  The 
theory  behind  the  technique  of  taking  visual  evoked  responses, 
is  that  the  chances  of  detecting  some  very  small  transient 
electromagnetic  field  may  be  enhanced  by  synchronizing  the 
sample  timing  so  that  50  to  a  100  of  these  synchronized  data 
sampling  windows  can  been  added  together,  and  then  averaged. 
Most  of  the  noise  signals  when  averaged  will  add  to  zero,  but 
if  there  is  a  good  signal  from  the  brain  which  is  hiding,  it 
will  grow  bigger,  since  it  will  be  synchronized. 

The  electroencephalograph  is  not  the  only  piece  of 
equipment  which  can  detect  visual  evoked  responses,  SQUID 
gradiometer  can  also  detect  them.  In  this  chapter,  we  will 
look  at  some  visual  evoked  responses  which  were  taken  with  the 
SQUID  gradiometer,  along  with  a  comparison  of  some  electro- 
encephalographic  visual  evoked  responses  verses  magnetoen- 
cephalographic  visual  evoked  responses.  For  further  research 
into  visual  evoked  responses  look  at  references  3,4,23,24,  and 


ment  was  taken  and  any  other  information  about  the 
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measurement  environment  may  be  there.  On  the  second 

line,  the  10  Hz  and  Notch  is  telling  us  that  the  10  Hz  low 
pass  filter  was  used,  along  with  the  60  Hz  notch  filter. 

Further  to  the  right  of  the  filter  information  is  the  value  of 

the  voltage  at  the  cursor  location.  Still  further  to  the  right 

is  the  vertical  scale  mode.  VLN  stands  for  the  vertical  scale 


is  in  linear  volts,  while  a  VLG  would  indicate  the  use  of  a 

logarithmic  scale.  On  the  third  line,  I  have  placed  the  date 

the  measurement  was  taken.  Sometimes,  especially  on  the  later 

plots,  I  entered  information  about  the  stimulus  on  this  line 

and  to  the  right  of  the  date.  Clear  to  the  far  right  of  the 

third  line  will  be  a  "T"  or  a  "C*.  For  all  of  the  measurements 

in  this  report,  a  "T"  will  be  to  the  far  right  of  the  third 

line  to  indicate  that  the  capture  control  setting  was  on 

transient.  If  I  wanted  to  take  measurements  with  the  capture 

control  setting  in  the  continuous  mode,  a  "C"  would  have 

replaced  the  "T".  On  the  fourth  line  you  will  find  the  cursor 

position  as  determined  by  the  horizontal  scale  units  selected 
for  the  measurement.  In  this  case,  the  cursor  position  is  in 
time,  that  is,  seconds.  Now  going  down  the  left  side  of  the 

plotted  box,  we  see  that  the  "A"  tells  us  that  the  input  of 

the  top  trace  is  was  on  channel  A,  while  the  "R"  tells  us  that 

the  format  of  that  trace  is  real,  versus  a  "P"  if  the  format 

was  phase,  or  a  "M"  if  the  format  was  magnitude,  or  "I"  if  the 

format  was  imaginary.  The  "TSU"  is  telling  us  what  type  of 

averaging,  while  the  "50"  is  telling  us  the  number  of  averages 

which  were  taken.  In  this  case,  the  "T*  stands  for  transient 

averaging,  while  the  "SU"  tells  us  the  sum  is  being  average. 

Other  possible  types  of  averaging  are  "PK"  for  peak,  "EX"  for 

expontential,  "DF"  for  difference,  and  "SW"  for  sweep. 

Continuing  down  the  left  side,  the  "AT"  tells  us  what  function 

was  performed.  In  this  case,  averaging  by  time  was  selected, 

but  other  types  of  function  could  have  been  selected.  For 

example,  a  "TF"  would  indicate  a  transfer  function  was 
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selected,  while  a  "RS"  would  indicate  that  a  RMS  spectrum 
function  was  selected.  Two  more  types  of  functions  are  "COH" 
for  coherence,  and  "it"  for  instantaneous  time.  If  we  were 
using  the  dual  trace  capability  of  this  FFT  analyzer,  simular 
information  on  the  input  channel  and  format  of  the  second 
trace  would  be  next  to  the  lower  left  side  of  the  box.  Now, 
going  across  the  bottom  of  the  box  from  left  to  right  we  see 
the  "2.0A"  and  the  "10.  A*  which  tells  us  the  maximum  setting 
for  channels  A  and  B,  respectively.  In  this  case,  the  maximum 
setting  for  channel  A  is  2  volts  ac,  while  the  maximum  setting 
for  channel  B  is  10  volts  ac.  A  "D"  instead  of  a  "A"  would 
have  indicated  the  values  were  dc  voltages  instead  of  ac 
voltages.  The  "B/16"  indicates  that  the  transient  trigger 
source  is  channel  B,  and  the  trigger  level  is  set  on  channel  B 
full  scale  divided  by  16.  Next,  the  "SEC"  indicates  that  the 
horizontal  scale  is  in  time,  that  is,  seconds.  The  horizontal 
scale  may  also  be  "HZ"  if  the  horizontal  scale  is  in 
frequency,  that  is,  hertz.  The  last  item  at  the  far  right  of 
the  bottom  line  is  the  horizontal  full  scale  value.  In  this 
case,  the  horizontal  full  scale  value  is  .4  second,  or  400 
milliseconds. 

Now  that  we  known  how  to  interpret  the  information  of 
the  plot,  lets  discuss  how  the  measurements  were  taken,  that 
is,  the  experimental  setup  for  the  human  magnetoencephalo- 
graphic  visual  evoked  response  experiment.  The  subject  sat  in 
a  wooden  chair  beneath  the  SQUID  gradiometer.  Pictures  of  the 
SQUID  gradiometer  will  be  shown  in  the  next  chapter.  The  SQUID 
gradiometer  mount  is  adjusted  so  that  the  SQUID  detection 


coil,  that  is,  the  bottom  of  the  dewar,  is  perpendicular  to, 
and  nearly  touching  the  place  on  the  subjects  head  at  which 
the  measurement  is  to  be  taken.  A  Grass  Model  PS22A  Photo¬ 
stimulator  was  used  to  trigger  a  bright  strobe  light  at  random 
intervals.  This  flash  was  used  as  the  stimulus.  The  signal 
used  to  fire  the  strobe  was  also  used  to  trigger  the  Nicolet 
660B  FFT  Analyzer.  For  the  visual  evoked  responses  measured 
from  humans,  the  FFT  analyzer  did  the  time  average  of  the  sum 
of  50  measurements,  each  of  a  400  millisecond  window  starting 
at  the  time  the  stimulus  was  triggered.  The  measurement  stored 
at  the  end  of  the  50  measurements  was  then  stored  on  an  8  inch 
disk  via  a  Intel  80/20  computer.  Later,  the  disk  files  were 
plotted  out  on  a  HP  digital  plotter. 


HUMAN  VER 

In  Figure  4-1,  is  an  example  of  a  human  visual  evoked 
response  taking  with  a  second-order  SQUID  gradiometer.  You 
will  find  this  plot,  and  other  human  visual  evoked  response 
plots  in  Appendix  C.  This  plot  is  a  typical  human  magnetoen- 
cephalographic  visual  evoked  response,  since  you  can  see  a  78 
millisecond  peak,  and  180  millisecond  peak.  These  latencies 
are  related  to  the  time  it  takes  the  signals  produced  in  the 
eye  to  be  conducted  via  the  optic  nerve,  pass  the  optic 
chiasm,  through  the  lateral  geniculate  nucleus,  and  finally  to 
the  primary  visual  cortex. 

As  you  look  thru  the  five  examples  of  human  magneto- 
encephalographic  visual  evoked  responses  which  are  contained 
in  Appendix  C,  note  that  all  five  different  human  visual 
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evoked  response  plots  had  an  initial  peak  around  80  milli¬ 
seconds.  This  component  of  the  magnetoencephalographic  visual 
evoked  response  seems  to  be  consistent. 


Figure  4-4  shows  us  the  EEG  profile  of  a  visual  evoked 
response . 


COMPONENTS  OF  THE  VISUAL  EVOKED  RESPONSE 
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FIGURE  4-4 
STANDARD  EEG  VER 


Note,  that  the  electroencephalographic  visual  evoked 
response  also  contains  a  peak  with  a  latency  off  approximately 
80  miliseconds,  along  with  a  peak  with  a  latency  off 
approximately  180  milliseconds. 


CHAPTER  V 


DOG  VISUAL  EVOKED  RESPONSE 


BACKGROUND 

During  the  late  part  of  September  and  early  part  of 
October  of  this  year,  we  got  the  opportunity  to  run  some 
magnetoencephalographic  visual  evoked  responses  experiments  on 
"Ricky",  the  laboratory  beagle  in  which  the  first  cortical 
brain  chip  was  implanted.  Appendix  D  has  a  number  of  plots  of 
"Ricky's"  magnetoencephalographic  visual  evoked  response.  We 
also  had  the  chance  to  take  some  magnetoencephalographic 
auditory  evoked  responses  from  Ricky.  The  plots  of  Ricky's 
auditory  evoked  responses,  along  with  some  discussion  of 
auditory  evoked  responses  are  contained  in  Appendix  E. 


EXPERIMENT  SETUP 

The  experiment  setup  for  the  Ricky's  magnetoencephalo¬ 
graphic  visual  evoked  response  measurements  were  similar  to 
the  setup  used  to  measure  the  human  magnetoencephalographic 
visual  evoked  responses.  For  example,  once  again  the  subject 
is  placed  beneath  the  SQUID  gradiometer,  and  the  SQUID  gradio- 
meter  mount  is  adjusted  so  that  the  SQUID  detection  coil,  that 
is,  the  bottom  of  the  dewar,  is  perpendicular  to,  and  nearly 
touching  the  place  on  the  subjects  head  at  which  the  measure¬ 
ment  is  to  be  taken.  Once  again,  a  Grass  Model  PS22A  Photo¬ 
stimulator  was  used  to  trigger  a  strobe  light  at  random 
intervals.  This  flash  was  used  as  the  stimulus.  The  signal 


used  to  fire  the  strobe  was  once  again  used  to  trigger  the  FFT 
analyzer.  Once  again,  the  Intel  80/20  computer  was  used  to 


store  the  average  of  the  50  measurements  on  a  8  inch  disk, 
which  was  used  later  to  recovery  the  data  so  that  it  could  be 
plotted  out.  The  major  difference  in  the  experiment  setup  for 
Ricky,  versus  the  experiment  setup  for  the  humans  is  that  the 
beagle,  Ricky,  was  lightly  anesthetized  to  prevent  him  from 
too  much  movement.  Pictures  of  the  experiment  setup  are  shown 
at  the  end  of  this  chapter. 

VER  COMPARISON 

In  Figures  5-1  and  5-2  are  two  of  the  magnetoencephalo- 
graphic  visual  evoked  response  which  were  taken  from  the 
laboratory  beagle,  Ricky.  While,  Figure  5-3  is  the  electroen- 
cephalographic  visual  evoked  response  taken  off  of  the  visual 
cortex  of  Ricky  by  the  AFIT  Semiconductor  Multielectrode 
Array,  which  was  implanted  on  his  brain.  We  expected  to  see  a 
magnetoencephalographic  visual  evoked  response  with  approxi¬ 
mately  the  same  latencies  as  the  latencies  measured  by  the 
AFIT  array.  Also,  as  we  move  the  position  at  which  the  SQUID 
gradiometer  is  measuring  the  visual  evoked  response,  we 
expected  to  see  a  phase  reversal.  It  turns  out  that  we  did  get 
approximately  the  same  latencies  as  did  the  AFIT  array.  Also, 
we  were  able  to  get  the  phase  reversal.  We  see  that  MEG. 022 
shown  in  Figure  5-1,  and  MEG. 045  shown  in  Figure  5-2  could 
possibly  be  a  phase  reversal,  but  the  best  example  of  a  phase 
reversal  is  from  two  of  the  plots  in  Appendix  D.  These  two 
plots  are  MEG. 057  and  MEG. 059. 
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FIGURE  5-4 
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FIGURE  5-8 


BARB  Me FAWN  WITH  STROBE  LIGHT 
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FIGURE  5-9 


ALVA  KARL  MONITORING  RICKY 


FIGURE  5-11 

DR.  HATSELL  POSITIONING  RICKY 
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CHAPTER  VI 


CONCLUSION 


DIPOLE  MODEL 

We  saw  that  the  predicted  values  in  both  the  magnetic 
field  and  the  surface  potential  had  differences,  especially  in 
absolute  value  in  relation  to  their  measured  values,  but  the 
wave  shape  of  the  measure  magnetic  field  or  surface  potential 
correlated  well  with  the  calculated  wave  shape,  except  for 
distortions  believed  to  be  due  to  such  factors  as  edge 
effects,  the  finite  volume,  and  the  lucite  mounting  block. 
Also,  the  conductivity  of  the  saline  solution  may  have  been  a 
factor  in  the  surface  potential  case. 


HUMAN  VER 

We  saw  some  excellent  plots  of  visual  evoked  responses 
in  humans,  long  with  some  discussion  into  VER  theory.  Also,  we 
saw  that  the  80  millisecond  component  of  the  visual  evoked 
response  taken  by  electroencephalographic  means,  also  shows  up 
in  the  magnetoencephalographic  visual  evoked  response. 

DOG  VER 


Once  again  we  saw  some  good  plots  of  visual  evoked 
responses  in  a  dog.  We  also  conclude  that  the  latency  of  the 
visual  evoked  response  as  measured  by  the  AFIT  implant  is 
approximately  the  same  as  the  latency  measured  by  the  SQUID 
gradiometer.  This  conclusion  would  tend  to  prove  that  the  AFIT 
array  seems  to  work. 


CHAPTER  VII 


RECOMMENDATIONS 


I  recommend  that  further  research  be  planned  especially 
with  a  researcher  looking  into  the  feasibility  of  simultaneous 
magnetoencephalographic  and  electroencephalographic  measure¬ 
ment.  Also,  more  research  be  done  in  the  area  of  determining 
means  by  which  the  SQUID  gradiometer  can  be  used  as  a 
diagnostic  tool.  For  example,  could  a  change  in  latencies  of 
various  evoked  responses  indicate  deterioration  of  the  visual 
system,  or  the  onset  of  Multiple  Sclerosis  (MS)  or  Alzheimer's 
disease? 
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fo  ISm  thesis  dipole  field  calculation  program 

20  REM  BFIELD  AND  VOLTAGE  (  X100  ) 

30  REM  BY  BOB  MURRAY 
40  REM 

50  LPRINT  CHR$ ( 12 ) 

60  INPUT  "ENTER  Z";Z 

70  LPRINT  TAB ( 15 ) ; " Z  =  ";Z 

80  LPRINT  TAB ( 15 ) ; " I  =  8  MILLIAMPS" 

90  LPRINT  TAB( 15 ) ; "L  =  .005  METERS" : LPRINT 

100  FOR  I=-10  TO  10 

110  FOR  J=-10  TO  10 

120  X=I/100 

130  Y=J/100 

140  T=4.D-12*X 

150  F=0 . 0 

160  B=( (X[2+Y[2+(.015+Z) [2)  [1.5) 

170  F=-T/B 
180  VF= ( 4 . 9D9 ) *F 

190  LPRINT  TAB (15) ;*B(";I;"," ; J; " ) =" ; F , TAB ( 15 ) ; " V( " ; I ; " , " ; J ; " ) =" ;VF 
200  NEXTJ 
210  NEXTI 

220  LPRINT  CHR$ ( 12 ) 

230  INPUT'ANOTHER  RUN";ANS$ 

240  IF( ANS$="Y" )  GOTO60 
250  END 


A-2 


1  -  8  MILL I AMPS 
L  «  .005  METERS 


B<-10 
B  ( -10 
B(  — 10 
B(-10 
BC-10 
B(  -10 
B<— 10 
B<-10 
B(-10 
B(  — 10 
B<— 10 
B<— 10 
B<  — 10 
B(  -10 
B(-10 
B<  -10 
B<  -10 
B<— 10 
B<— 10 
B<— 10 
B(— 10 
B<  — 9  , 
B<  -9  , 
B<— 9  , 
B<— 9  , 
B<— 9  , 
B<— 9  , 
B(— 9  , 
B(— 9  , 
B(-9  , 
B(-9  , 
B<— 9  , 
B<  -9  , 
B<— 9  , 
B<— 9  , 
B<— 9  , 
B<— 9  , 
B(— 9  , 
B(— 9  , 
B(-9  , 
B(— 9  , 
B<— 9  , 
B(-8  , 
B<-8  , 
B<-8  , 
B<-8  , 
B(-B  , 
B<-8  , 
B(-8  , 
B<-8  , 
B<-8  , 
B<-8  , 
B<  -8  , 
B<-8  , 
B<-8  , 
B(— 8  , 
B<-8  , 
B(— 8  , 
B(-8  , 
B(-8  , 
B(-B  , 


,-10  )-  1.28051E-10 

V(-10 

,-10  ) 

=  .627451 

• -9  )■  1.47246E-10 

V(-10 

, -9  )  = 

.721507 

, -8  )=  1.68875E-10 

V(-10 

,  -8  >  = 

. 827488 

, -7  >=  1 . 92761E-10 

V(-10 

,  -7  )  = 

.944528 

, -6  )=  2. 1B409E-10 

V(-10 

,-6  )  = 

1 . 0702 

, -5  )=  2. 44902E-10 

V<-10 

,-5  )  = 

1 . 20002 

, -4  )=  2.70832E-10 

V(-10 

,  -4  >  = 

1 . 32708 

, -3  )=  2. 94338E-10 

V<-10 

,  -3  )* 

1 . 44226 

, -2  )=  3.13293E-10 

V(-10 

,  -2  )  = 

1.53514 

,-l  )=  3. 25665E-10 

V<-10 

,-l  >  = 

1.59576 

,  0  )=  3.29972E-10 

Vl-10 

,  0  )  = 

1.61686 

,  1  )=  3.25665E-10 

V<-10 

,  1  )  = 

1.59576 

,  2  )=  3. 1 3293E-10 

V(-10 

,  2  >  = 

1.53514 

,  3  )=  2.94338E-10 

V<-10 

,  3  )  = 

1.44226 

,  4  )  =  2. 70832E-10 

V(-10 

,  4  )  = 

1 . 32708 

,  5  )  =  2.44902E-10 

V(-10 

,  5  >  = 

1 . 20002 

,  6  )=  2.18409E-10 

V(-10 

,  6  )  = 

1 . 0702 

,  7  )=  1 . 92761E-10 

V(-10 

,  7  >  = 

. 944528 

,  8  )=  1 . 68875E-10 

V(-10 

,  8  )  = 

. 827488 

,  9  >=  1 . 47246E-10 

V(-10 

,  9  )  = 

.721507 

,  10  >  =  1 . 2805 IE- 10 

V<-10 

,  10  ) 

-  .627451 

-10  )=  1 . 32522E-10 

V(-9 

,-10  >= 

. 649356 

-9  )=  1.5459E-10 

V(-9 

,  -9  >  = 

. 757492 

-8  )=  1 . 80085E— 10 

V(-9 

,  -8  )  = 

.882418 

-7  )=  2. 09008E-10 

V(-9 

,  -7  )  = 

1.02414 

-6  )=  2.40957E-1B 

V<-9 

,-6  >- 

1.18069 

-5  )=  2. 74926E— 10 

V(-9 

, -5  )* 

1.34714 

-4  )=  3. 09129E— 10 

V(-9 

,  —4  >  = 

1.51473 

-3  >«  3.40955E-10 

V(— 9 

,  -3  )  = 

1.67068 

-2  )=  3.67192E— 10 

V(-9 

,  -2  >  = 

1 . 79924 

-1  )=  3.84594E-10 

V<-9 

,  -1  )  = 

1.88451 

0  )=  3.90703E-10 

V(— 9 

,  0  )* 

1.91444 

1  )=  3. 84594E— 10 

V(  -9 

,  1  )  = 

1.88451 

2  )=  3. 67192E— 10 

V(— 9 

,  2  )  = 

1 . 79924 

3  )»  3. 40955E-10 

V(-9 

,  3  )  = 

1 . 67068 

4  )*  3.09129E-10 

V(-9 

,  4  )  = 

1.51473 

5  )=  2. 74926E-10 

V(-9 

,  5  >  = 

1.34714 

6  )=  2. 40957E— 10 

V(-9 

,  6  )  = 

1.18069 

7  )=  2. 09008E-10 

V(-9 

,  7  )  = 

1.02414 

8  >*  1 . 80085E— 10 

V(-9 

,  8  )« 

.882418 

9  )=  1 . 5459E-10 

V<-9 

,  9  >  = 

.757492 

10  )  =  1 . 32522E-10 

V(-9 

,  10  )  = 

. 649356 

-10  )=  1 . 35 IE— 10 

V(-8 

,-10  )  = 

.66199 

-9  >=■  1.60076E-10 

V(  -8 

,-9  )  = 

. 784371 

-8  >=  1 . 89732E— 10 

V(-8 

,  -8  >« 

. 929687 

-7  )-  2. 24407E-10 

V<-8 

,  -7  >« 

1.09959 

-6  >-  2. 63977E-10 

V(-8 

, -6  >- 

1 . 29349 

-5  )=  3. 07509E-10 

V(-8 

, -5  >« 

1 . 50679 

-4  )=  3.52866E-10 

V(-8 

,-4  )- 

1 . 72904 

-3  )«  3. 96457E-10 

V(-8 

, -3  >  = 

1 . 94264 

-2  >«  4. 33408E-10 

V(-8 

,  -2  >- 
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-1  )»  4. 58428E-10 

V(-8 

,-i  >■ 

2.2463 
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V(  -8 

,  0  )■ 
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1  )>  4.58428E-10 

V(-8 

,  1  >« 
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2  >-  4. 3340BE-10 

V(-8 

,  2  )- 
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3  )-  3. 96457E-10 

V<-8 

,  3  >« 
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4  >-  3.52866E-10 

V(-8 

,  4  >“ 

1 . 72904 

5  >*  3. 07509E-10 

V(-B 

,  5  )- 

1.50679 

6  )«  2. 63977E-10 

V(  -8 

,  6  >* 

1 . 29349 

7  )-  2. 24407E-10 

V<-8 

,  7  >« 
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8  )=  1 . 89732E-10  , 

V(-8 

,  .8  >■ 

.929687  _ 

.y 


B(  — 8 

9 

)  = 

1 . A0076E-10 

VC-  R 

.  9 

)  = 

. 78437 t  - 

BC-C 

10 

>  = 

1 . 35 IE- 10 

VC-8 

.  10 

>=  .66199 

B(-7 

-10 

>  = 

1 . 34933E-10 

V(-7 

,-10 

)  = 

.661169 

B<-7 

-9 

)  = 

1.62562E-10 

VC  -7 

,  “9 

)  = 

. 796552 

BC-7 

-8 

)® 

1.96356E-10 

VC-7 

,-8 

)« 

.962143 

B(-7 

-7 

)  = 

2. 37212E-10 

VC— 7 

,-7 

)« 

1.16234 

B<-7 

-6 

)  = 

2.85594E-10 

VC-7 

*  —6 

)  = 

1.39941 

B<-7 

-5 

)« 

3.40983E-10 

VC-7 

,-5 

)  = 

1 . 67082 

B<-7 

-4 

)  = 

4.01 124E— 10 

VC-7 

,-4 

)  = 

1.96551 

B<-7 

-3 

)  = 

4.61286E-10 

VC-7 

,-3 

>*» 

2.2603 

B<-7 

-2 

)  = 

5. 14 123E-10 

VC-7 

,-2 

>  = 

2.5192 

B<-7 

-1 

)  = 

5.50873E-10 

VC-7 

,-l 

)  = 

2.69928 

B<— 7 

0 

)  = 

5.64106E-10 

VC-7 

,  0 

)  = 

2.76412 

B<-7 

1 

>* 

5.50873E-10 

VC-7 

,  1 

)  = 

2.69928 

B(-7 

2 

>  = 

5. 14123E-10 

VC-7 

,  2 

)s 

2.5192 

B<-7 

3 

>  = 

4.612B6E-10 

VC  -7 

,  3 

)*= 

2.2603 

B(-7 

4 

)* 

4.01 124E-10 

VC  -7 

.  4 

)= 

1.96551 

B(— 7 

5 

)* 

3. 40983E-10 

VC  -7 

,  5 

)= 

1 . 67082 

B(-7 

6 

>  = 

2.85594E-10 

VC  -7 

*  6 

)* 

1.39941 

B<-7 

7 

)  = 

2. 37212E-10 

VC-7 

,  7 

)= 

1 . 16234 

BC-7 

8 

)  = 

1 . 96356E-10 

VC-7 

,  8 

)= 

.962145 

B<  — 7 

9 

)* 

1 . 62562E-10 

VC  -7 

,  9 

)= 

. 796552 

B<-7 

10 

)«  1 . 34933E-10 

VC— 7 

,  10 

)= 

.661169 

B<-6 

-10 

)=  1 . 31043E-10 

VC-6 

,-10 

)- 

.642123 

B<  — 6 

-9 

)  = 

1.60638E-10 

VC -6 

,  —9 

)* 

.787125 

B(— 6 

-8 

)  = 

1 . 97983E-10 

VC-6 

,  -8 

)= 

. 9701 16 

B(— 6 

-7 

>  = 

2. 44794E— 10 

V<-6 

,  -7 

)= 

1.19949 

B(  — 6 

-6 

)  = 

3.02563E-10 

VC— 6 

,  —6 

)= 

1.48256 

B(— 6 

-5 

)  * 

3. 71807E-10 

VC-6 

,-5 

)= 

1.82185 

FC-6 

)  = 

4. 50778E— 10 

VC— 6 

,  -4 

)= 

2.20881 

B<  -6 

-3 

)  = 

5 . 33783E— 10 

V(— 6 

,-3 

)  = 

2.61554 

BC-6 

-2 

)* 

6. 10058E-10 

VC-6 

i-2 

)  = 

2.98928 

BC-6 

-1 

)  = 

6.65009E-10 

VC— 6 

,-l 

)= 

3.25853 

BC-6 

0 

)  = 

6.85185E-10 

VC— 6 

,  0 

>“ 

3.33741  i 

B<— 6 

1 

)* 

6. 65009E-10 

V(— 6 

*  1 

)  = 

3.25833 

B(— 6 

2 

)  = 

6. 10058E— 10 

V(— 6 

,  2 

)= 

2.98928 

BC— 6 

3 

)« 

5.33783E-10 

VC -6 

,  3 

>= 

2.61534 

B<— 6 

4 

)  = 

4. 50778E— 10 

VC— 6 

,  4 

)  = 

2.20881 

BC-6 

5 

)  = 

3.71807E-10 

VC— 6 

,  5 

)= 

1.82185 

B<— 6 

6 

)» 

3.02563E-10 

VC -6 

t  6 

)= 

1.48256 

BC-6 

7 

)  = 

2. 44794E-10 

VC -6 

,  7 

)= 

1.19949 

B<— 6 

8 

)* 

1.97983E-10 

VC-6 

,  8 

>* 

.970116 

BC-6 

9 

)» 

1 . 60638E-10 

VC -6 

,  9 

)» 

.787123 

B<-6 

10 

)»  1 .31045E-10 

VC-6 

,  10 

)>  .642123 

B<— 5 

-10 

)=  1.22431E-10 

VC -5 

,-10 

)*  .60001 

B(-5 

-9 

)  = 

1.52737E-10 

VC -5 

,  -9 

)  = 

. 748409 

B(— 3 

-8 

1 . 92193E-10 

VC -5 

,-8 

)» 

.941746 

B<-5 

-7 

)  = 

2.4356E-10 

VC-3 

,-7 

)  = 

1.19344 

BC-5 

-6 

>  = 

3.09839E-10 

VC -5 

v  —6 

)  = 

1.51821 

B(-5 

-5 

)* 

3. 9348E-10 

VC-5 

»-5 

)  = 

1.92803 

B<-3 

-4 

)« 

4.94502E-10 

VC-5 

,  -4 

>  = 

2.42306 

B(-3 

-3 

)* 

6.0728E-10 

VC-5 

,-3 

>=* 

2.97367 

B(-5 

-2 

)  = 

7.17036E-10 

VC-5 

,-2 

>® 

3.51348 

BC-5 

-1 

)  = 

7.99B5E-10 

VC-5 

,-l 

>=* 

3.91926 

B<-5 

0 

)  = 

8.31059E-10 

VC-5 

,  0 

)  = 

4.07219 

B(-3 

1 

>* 

7.9985E-10 

VC -5 

,  1 

)  = 

3.91926 

B<-5 

2 

>* 

7 . 17036E-10 

VC -3 

,  2 

)  = 

3.51348 

B<-3 

3 

)* 

6. 0728E-10 

VC-5 

,  3 

>» 

2.97367 

BC-5 

4 

>* 

4.94502E-10 

VC-5 

,  4 

)- 

2.42306 

BC-5 

5 

)* 

3. 9348E-10 

VC -5 

»  5 

)  = 

1 . 92803 

BC-5 

6 

)* 

3.09839E-10 

VC-5 

t  6 

)  = 

1.31821 

BC-5 

7 

)* 

2. 4356E-10 

VC -3 

,  7 

)» 

1.19344 

B(-5 

8 

)- 

1.92193E-10 

VC-5 

,  8 

)=» 

.941746 

B<-5 

9 

>* 

1 . 52737E-10 

VC-5 

,  9 

)  = 

. 748409 

B(-3 

10 

)=  1 . 22431E-10 

VC-5 

,  10 

>»  .60001 

BCr4 

.-10  J.1 

*  1 . 08333E-10 

..VC -4. 

-u-10-1.* 

=_  .530831 

A-4 
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BC-4 

-9 

>  = 

1 . 37391E-10  “ 

V<-4 

-9 

)- 

.673214 

B<-4 

-8 

)* 

1.76433E-10 

VC  -4 

-8 

)  = 

.864522 

BC— 4 

-7 

>  = 

2.29214E-10 

VC-* 

-7 

><* 

1.12313 

B(  -4 

-6 

)  = 

3.00518E-10 

VC -4 

-6 

)  = 

1.47254 

B<-4 

-5 

)  = 

3.95602E-10 

VC-4 

-5 

)  = 

1 . 93843 

B<-4 

-4 

>  = 

5 . 1807E-10 

VC -4 

-4 

)* 

2.53854 

B(  -4 

-3 

>=» 

6. 64847E-10 

V(  -4 

-3 

>  = 

3.25773 

B<— 4 

-2 

)BS 

8. 18217E-10 

VC  -4 

-2 

)  = 

4 . 00926 

BC-4 

-1 

)  = 

9. 41073E-10  . 

VC— 4 

-1 

)  = 

4.61127 

BC— 4 

0 

>  = 

9. 89019E-10 

VC -4 

0 

)  = 

4.84619 

BC-4 

1 

)  = 

9. 41075E-10 

VC-4 

1 

)  = 

4.61127 

B(-4 

2 

)* 

8.18217 E-10 

VC-4 

2 

)=■ 

4.00926 

B(  -4 

3 

)  = 

6. 64847E-10 

VC-4 

3 

)  = 

3.25773 

B(  -4 

4 

)  = 

5. 1807E-10 

VC-4 

4 

)  = 

2.53834 

B<-4 

5 

)  = 

3.95602E-10 

VC-4 

5 

)  = 

1 . 93843 

B<-4 

6 

)  = 

3. 005 18E-10 

VC-4 

6 

)  = 

1 . 47254 

B<  — 4 

7 

>  = 

2.29214E-10 

VC -4 

7 

)  = 

1.12313 

B(— 4 

8 

)  = 

1 . 76433E— 10 

VC-4 

8 

>  = 

. 864522 

BC— 4 

9 

)  = 

1 . 37391E-10 

VC  -4 

9 

>  = 

.673214 

BC— 4 

10 

)  = 

1 . 08333E-10 

VC -4 

10 

>  = 

.530831 

B<-3 

-10 

)=  8. 83014E-1 1 

VC-3 

-10 

>=  .432677 

BC— 3 

-9 

)  = 

1. 13652E-10 

VC-3 

-9 

)  = 

. 556894 

B<-3 

-8 

)  = 

1 . 48671E-10 

VC -3 

-8 

>  = 

. 72849 

B<-3 

-7 

)  = 

1.97694E-10 

VC -3 

-7 

)= 

. 968701 

B<-3 

-6 

)  = 

2. 66892E-10 

VC  -3 

-6 

)  = 

1.30777 

BC-3 

-5 

)  = 

3. 64368E-10 

VC-3 

-5 

>  = 

1 . 7834 

B<— 3 

-4 

)= 

4.98633E-10 

VC-3 

-4 

)  = 

2. 44331 

BC-3 

-3 

)* 

6.72663E-10 

VC-3 

-3 

)  = 

3.29603 

B<-3 

-2 

)  = 

8. 70279E-10 

VC-3 

-2 

)  = 

4.26437 

B<-3 

-1 

)s 

1.04072E-09 

VC-3 

-1 

)  = 

5.09952 

BC— 3 

0 

)  = 

1 . 1 1027E-09 

VC-3 

0 

)  = 

5.44033 

B(-3 

1 

)  = 

1.04072E-09 

VC-3 

1 

><= 

5.09932 

B(-3 

2 

)= 

8.70279E-10 

VC-3 

2 

)  = 

4.26437 

B<-3 

3 

)a 

6.72663E-10 

VC-3 

3 

)» 

3.29603 

B(— 3 

4 

)=> 

4. 98635E-10 

VC -3 

4 

>® 

2.44331 

BC— 3 

5 

)* 

3. 64368E-10 

VC-3 

5 

)  = 
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B< 

8 

,  ~4 

)=-3.52866E-10 

v< 

8 

,"4 

>=-1.72904 

B< 

8 

,  “3 

)=-3.96457E-10 

V( 

8 

,  -3 

>=-1.94264 

B( 

8 

,  ”2 

)=-4. 33408E-10 

V( 

8 

,-2 

>=-2. 1237 

B( 

8 

,-l 

)— 4.58428E-10 

V( 

8 

,-l 

>=-2.2463 

B( 

8 

»  0 

)  *-4 . 67307E- 1 0 

V( 

8 

,  0 

>=-2.2898 

B< 

8 

,  1 

)=-4. 58428E-10 

v< 

8 

,  1 

>=-2.2463 

B< 

8 

,.  2. 

)=-4.33408E.“J0_  _  _ 

8 

-t.2. 

) =-2 . 1237 _ 

B( 

8  ,  3  )=-3.96437E-10 

V( 

8  , 

3  >=-1.94264' 

B( 

8  ,  4  >=-3.52866E-10 

V( 

8  , 

4  >=-1.72904 

B( 

8  ,  5  >=-3.07509E-10 

V< 

8  , 

5  >=-1.50679 

B( 

8  ,  6  >=-2. 63977E-10 

V( 

8  , 

6  >=-1.29349 

B( 

8  ,  7  )=-2. 24407E— 10 

V( 

8  , 

7  >=-1.09959 

B< 

8  ,  8  )  =- 1 . 89732E- 10 

V< 

8  , 

8  >=-.929687 

B( 

8  ,  9  )=-l . 60076E-10 

V( 

8  , 

9  >=-.784371 

B( 

8  ,  10  >=— 1.351E-10 

V( 

8  , 

10  >=-.66199 

B( 

9  ,-10  )=-1.32522E-10 

V( 

9  , 

-10  >=-.649356 

B( 

9  i -9  >=-l .5459E-10 

V< 

9  , 

-9  >=-.757492 

B( 

9  , -8  )=-l . 80085E-10 

V( 

9  , 

-8  >=-.882418 

B( 

9  , -7  ) =-2. 09008E-10 

V( 

9  , 

-7  >=-1.02414 

B( 

9  , -6  )=-2. 40957E-10 

V( 

9  , 

-6  >=-1.18069 

B( 

9  , -5  >=-2.74926 E-10 

V( 

9  , 

-5  >=-1.34714 

B< 

9  , -4  )=-3. 09129E-10 

V< 

9  , 

-4  >=-1.51473 

B< 

9  , -3  )=-3. 40955E-10 

V( 

9  , 

-3  >=-1.67068 

B( 

9  , -2  )=-3. 67192E-10 

V< 

9  , 

-2  >=-1.79924 

B( 

9  ,-l  ) =-3. 84594E-10 

9  , 

-1  >=-1.88451 

B( 

9  ,  0  ) =-3. 90703E-10 

V( 

9  , 

0  >=-1.91444 

B( 

9  ,  1  )=-3. 84594E-10 

V< 

9  , 

1  >=-1.88451 

B( 

9  ,  2  >=-3. 67192E-10 

V< 

9  , 

2  >=-1.79924 

B< 

9  ,  3  >— 3.40955E-10 

V( 

9  , 

3  >=-1.67068 

B( 

9  ,  4  )=-3.09129E-10 

V< 

9  , 

4  >=-1.51473 

B( 

9  ,  5  >=-2. 74926E— 10 

V< 

9  , 

5  >=-1.34714 

B( 

9  ,  6  ) =-2. 40957E— 10 

V< 

9  , 

6  >=-1.18069 

B( 

9  ,  7  )=-2. 09008E— 10 

V< 

9  , 

7  >=-1.02414 

B( 

9  ,  8  ) =- 1 . 80085E— 1 0 

V< 

9  , 

8  >=-.882418 

B( 

9  ,  9  )=-l . 5459E-10 

V( 

9  , 

9  >=-.757492 

B( 

9  ,  10  )  =- 1 . 32S22E- 1 0 

VC 

9  , 

10  >=-.649356 

B( 

10  ,-10  )=-l . 28051E-10 

vc 

10 

,-10  >=-.627451 

B( 

10  , -9  >=-l . 47246E— 10 

VC 

10 

,-9  >=-.721307 

B( 

10  , -8  )*-l . 68875E-10 

10 

, -8  >=-.827488 

B( 

10  , -7  )»-l . 92761E— 10 

V( 

10 

, -7  >=-.944528 

B( 

10  , -6  )=-2. 18409E-10 

v< 

10 

, -6  >=-1.0702 

B< 

10  , -5  >=-2.44902E-10 

V( 

10 

, -5  >=-1.20002 

B< 

10  , -4  )=-2. 70832E— 10 

v< 

10 

» -4  >=-1.32708 

B< 

10  , -3  >=-2. 94338E-10 

V( 

10 

, -3  >=-1.44226 

B< 

10  , -2  )=-3. 13293E-10 

V( 

10 

, -2  >=-1.53314 

B( 

10  ,-l  )=-3. 25665E-10 

V( 

10 

,-l  >*-1.59376 

B( 

10  ,  0  )=-3. 29972E-10 

v< 

10 

,  0  >=-1.61686 

B( 

10  ,  1  )=— 3. 2S665E-10 

V( 

10 

,  1  >=-1.59576 

B( 

10  ,  2  )=-3. 13293E-10 

v< 

10 

,  2  >*-1.53514 

B  < 

10  ,  3  >— 2.94338E-10 

v< 

10 

,  3  >=-1.44226 

B( 

10  ,  4  ) =-2 . 70832E- 1 0 

V( 

10 

,  4  >=-1.32708 

B( 

10  ,  5  )=-2. 44902E-10 

v< 

10 

,  5  >=-1.20002 

B( 

10  ,  6  )*-2. 18409E-10 

V( 

10 

,  6  >=-1.0702 

B( 

10  ,  7  >=-l . 92761E-10 

v< 

10 

,  7  >=-.944528 

B( 

10  ,  8  >«- 1.68875E-10 

V( 

10 

,  8  >=-.827488 

B( 

10  ,  9  )=-l . 47246E-10 

V( 

10 

,  9  >=-.721307 

B( 

10  ,  10  )=-1.28051E-10 

v< 

10 

,  10  >=-.627431 

TABLE  A-l 


Comparison  of  Predicted  Versus  Measured  Magnetic  Field 


Coordinates 

Predicted 

Measured 

Facte 

(x,y) 

mV 

mV 

-10,-10 

0.627 

.042 

14.9 

-9,-9 

0.757 

.049 

15.4 

-8,-8 

0.929 

.063 

14.7 

-7,-7 

1.162 

.081 

14.3 

-6,-6 

1.483 

.126 

11.8 

-5,-5 

1.928 

.234 

8.2 

-4,-4 

2.539 

.371 

6.8 

-3,-3 

3.296 

.61 

5.4 

-2,-2 

3.881 

.788 

4.9 

-1,-1 

3.154 

.72 

4.4 

o 

o 

0.000 

.4 

1,1 

3.154- 

.731 

4.3 

2,2 

3.881- 

.811 

4.8 

3,3 

3.296- 

.633 

5.2 

4,4 

2.539- 

.387 

6.6 

5,5 

1.928- 

.251 

7.7 

6,6 

1.483- 

.147 

10.1 

7,7 

1.162- 

.088 

13.2 

8,8 

0.929- 

.067 

13.9 

9,9 

0.757 

.051 

14.8 

10,10 

0.627 

.043 

14.6 

10  REM 

20  REM  PROGRAM  PHIB/BAS 

30  REM  CALCULATES  SURFACE  POTENTIAL  FOR  A  DIPOLE 

40  REM  SUBMERGED  IN  A  .464%  NACL  SOLUTION 

50  REM  1=8  MILLIAMPS 

60  REM  LENGTH  OF  DIPOLE  =  .005  METERS 

70  REM  CAPT.  ROBERT  D.  MURRAY 

80  REM 

90  P=. 008*. 005 

100  PI=3. 14159 

110  D= . 0 1 5 

120  SIGMA= .7624 

130  FOR  X=- . 1  TO  .1  STEP  .01 

140  FOR  Y=- . 1  TO  .1  STEP  .01 

150  Rl=( (X*X+Y*Y+D*D) [1.5) 

160  R2=( (X*X+Y*Y+D*D) [1.5) 

170  K=P/(2*PI*SIGMA) 

180  PHI  =  ( (K*Y)*( (1/R1)+{1/R2) ) ) 

190  LPRINT"PHI ( " ; INT(100*( . 005+X) )?","? INT(100*( .005+Y) ) ; " ) 
200  NEXTY 
210  NEXTX 


PHI  <-10  ,-10  )  -  -5.80624E-04 
PHI <-10  i -9  )  ■  -6. 05905E-04 
PHI <-10  , -8  )  =  -6. 23268E-04 
PHI <-10  ,-7  )  =  -6.28467E-04 
PHI <-10  , -6  )  =  -6.16426E-04 
PHI <-10  , -5  )  =  -5. 81716E-04 
PHI <-10  , -4  )  =  -5.195E-04 
PHI <-10  , -3  )  =  -4.26971E-04 
PHI <-10  »— 2  )  =  -3. 04975E-04 
PHI <-10  ,-l  )  »  -1 . 59182E-04 
PHI <-10  ,  0  >  =  -2. 40687E— 10 
PHI <-10  ,  1  >  *  1 . 59181E-04 
PHI <-10  ,  2  >  =  3. 04973E-04 
PHI <-10  ,  3  >  =  4. 26971E-04 

PHI <-10  ,  4  )  =  5.195E-04 

PHI <-10  ,  5  )  =  5.8171 6E-04 

PHI <-10  ,  6  >  =  6. 16426E-04 

PHI <-10  ,  7  )  =  6.28467E-04 

PHI <-10  ,  8  >  =  6.23269E-04 

PHI <-10  ,  9  >  x  6.05905E-04 
PHI <-10  ,  10  >  =  5 . 80623E-04 

PHI <-9  ,-10  >  =  -6. 73227E-04 
PHI <-9  , -9  >  =  -7. 14023E-04 
PH I <  — 9  , -8  )  =  -7. 47712E-04 
PH I <  — 9  , -7  )  =  -7.68636E-04 
PH I < -9  , -6  )  =  -7.69469E-04 
PH I < -9  , -5  )  =  — 7. 41405E-04 
PH I <  —9  , -4  )  =  -6.75604E-04 
PH I < -9  ,-3  >  =  —5 . 65457E— 04 
PHI <-9  , -2  )  »  -4 . 09836E— 04 
PH I < -9  ,-l  )  *  -2. 1596E-04 
PHI <— 9  ,  0  )  *  -3. 27621E-10 
PHI <-9  ,  1  )  =  2. 15959E-04 

PH I < -9  ,  2  )  =  4. 09835E-04 

PH  I  <  -9  ,  3  )  »  5 . 65456E-04 

PHI <-9  ,  4  >  =  6.75603E-04 

PHI <-9  ,  5  )  =  7. 41405E-04 

PHI <-9  ,  6  )  o  7 . 69469E-04 
PHI <-9  ,  7  >  -  7.68656E-04 

PH I < -9  ,  8  )  -  7. 47712E-04 

PHI <-9  ,  9  )  *  7. 14023E-04 

PH I < -9  ,  10  )  =  6 . 73227E-04 

PHI <— 8  ,-10  )  =  -7 . 79085E-04 
PHI <-8  , -9  >  *  -8.4 1 176E-04 
PHI <-8  ,-B  )  *  -8. 98774E-04 
PHI <-8  ,-7  )  ■  -9. 44853E-04 
PHI <-B  , -6  )  =  -9. 69134E-04 
PHI <-8  , -5  )  =  -9. 57962E-04 
PHI <-8  , -4  )  ■  -8. 95537E-04 
PHI <-8  , -3  >  =  -7. 67519E-04 
PHI <-8  , -2  )  =  -5 . 67268E-04 
PH I <  — 8  ,-l  )  =  -3. 02824E-04 
PH I < -8  ,  0  )  ■  -4. 61497E— 10 
PHI <-8  ,  1  )  ■  3. 02823E-04 
PHI <-8  ,  2  )  ■  5 . 67267E-04 
PHI<-8  ,  3  )  «  7 . 675 1 8E-04 

PHI<-8  ,  4  )  «  8. 95537E-04 

PHI <-B  ,  5  )  «  9. 57962E-04 

PH I < -8  ,  6  )  ■  9. 69134E-04 

PH I < -8  ,  7  )  -  9. 44853E-04 

PHI <-8  ,  8  )  =  8. 98775E-04 

PHI < -8  ,  9  )  «  8. 41 176E-04 

PHI <-B  ,  10  >  -  7.79086E-04 

PH I < -7  , - 1 0_ )  »  -8 .9781 E-04 _ 


p,7  V  >  V  T*  *7  V  •  V  V 

•;•  'V  V\v*>  V-  ■» 

C-> 

E 

PHI <-7 

i -9  ) 

= 

-9. 83272E-04 

UN- 

Ki  j:* 

/*,%  w 

PHI < -7 

.-8  ) 

s 

-1.07983E-03 

PHI < -7 

i-7  ) 

= 

-1. 16466E-03 

r.‘l 

PHI (-7 

i-6  ) 

= 

-1.22951E-03 

H 

PHI (-7 

.-5  > 

as. 

-1 • 25412E-03 

. 

PHI (-7 

,-4  ) 

= 

-1. 21 129E-03 

PHI (-7 

.-3  ) 

= 

-1. 07131E-03 

::::  C 

•  / 

PHI (-7 

i-2  ) 

= 

-8. 13316E-04 

PHI (-7 

i-l  ) 

= 

-4. 42179E-04 

PHI (-7 

i  0  ) 

= 

-6. 78276E-10 

PHI (-7 

i  1  > 

= 

4. 42178E-04 

■ 

PHI (-7 

i  2  ) 

= 

8. 13315E-04 

w. 

PHI <-7 

i  3  ) 

= 

1. 07131E-03 

*.-. 

PHI (-7 

i  4  ) 

= 

1.211 29E-03 

>;■ 

PHI (-7 

i  5  ) 

sx 

1.25412E-03 

•*. 

PHI (-7 

i  6  ) 

= 

1 . 2295  IE-03 

PHI (-7 

i  7  ) 

= 

1. 16466E-03 

PHI <-7 

i  8  ) 

= 

1 . 07983E-03 

■ 

PHI <-7 

i  9  ) 

= 

9. 88272E-04 

PHI <-7 

,  10 

)  = 

8.9781E-04 

:< 

PHI<-6 

,-10 

)  = 

-1 . 02738E-03 

PHI (-6 

i-9  > 

= 

-1. 1542E-03 

PHI (-6 

i-8  ) 

s 

-1.29218E-03 

S; 

PHI<-6 

i-7  ) 

= 

-1 . 43442E-03 

>; 

PHI (-6 

i-6  ) 

= 

-1.56615E-03 

fell] 

PHI <-6 

i-5  ) 

= 

-1.65999E-03 

Pf 

PHI <-6 

i-4  ) 

= 

-1 . 67181E-03 

SN, 

PHI <-6 

i-3  ) 

= 

-1 . 54257E-03 

£■•’ 

PHI (-6 

,-2  > 

= 

-1.21623E-03 

fcV 

PHI <-6 

i-l  ) 

= 

-6. 79155E-04 

p. 

PHI (-6 

i  0  > 

= 

-1 . 05196E-09 

PHI (-6 

i  1  > 

= 

6. 79153E-04 

1  **  c 

PHI (-6 

i  2  > 

3 

1.21623E-03 

PHI (-6 

i  3  ) 

3 

1.54237E-03 

> 

PHI <-6 

i  4  ) 

3 

1.67181E-03 

5:- 

PHI <-6 

i  5  ) 

3 

1 . 65999E-03 

PHI <-6 

i  6  ) 

= 

1 . 56615E-03 

«/•  ’ 

PH I < -6 

i  7  ) 

a 

1 . 43442E-03 

•si 

PHI (-6 

i  8  ) 

3 

1.29218E-03 

■ 

PHI (-6 

i  9  ) 

3 

1. 1542E-03 

i 

PHI (-6 

i  10 

)  a 

1.02738E-03 

>: . 

PHI (-5 

,-10 

)  * 

-1.16343E-03 

& 

PHI (-5 

i-9  ) 

m 

-1.33453E-03 

PHI (-5 

,-8  ) 

3 

-1 . 53274E-03 

■‘w 

PHI (-5 

i-7  ) 

* 

-1.75576E-03 

‘•V  ■ 

PHI (-5 

i-6  ) 

s 

-1.99199E-03  4 

k  ., 

PHI (-3 

i-5  ) 

3 

-2. 21089E-03 

1  ! 

PHI (-3 

i-4  ) 

3 

-2.3486E-03 

PHI (-5 

i-3  ) 

3 

-2. 29555E-03 

>: 

PHI  C-3 

i-2  ) 

3 

-1 . 91 198E-03 

$ 

PHI (-5 

i-l  > 

3 

-1.11224E-03 

$> 

PHI (-5 

i  0  > 

3 

-1 . 74943E-09 

s* 

PHI (-5 

i  1  > 

3 

1.  U224E-03 

PHI <-3 

i  2  ) 

8 

1 . 91 198E-03 

m . 

PHI (-5 

i  3  ) 

8 

2. 29555E-03 

PHI <-5 

i  4  ) 

= 

2. 3486E-03 

PHI <-5 

i  5  > 

a 

2.2109E-03 

,*  * 

^  o 

u  ^ 

PHI <-5 

i  6  ) 

3 

1 . 99199E-03 

PHI (-5 

i  7  ) 

8 

1 . 73576E-03 

PHI (-3 

i  8  ) 

3 

1 . 53274E-03 

a 

PHI <-3 

i  9  ) 

S 

1.33453E-03 

1 

PHI <-5 

i  10 

>  = 

1. 16343E-03 

PHI < -4 

,-10 

)  - 

-1.29875E-03 

:*:  ;v 

PHI <-4 

i-9  ) 

3 

-1.3201  IE-03 

>.  V 

PHI (-4 

i-8  ) 

3 

-I.79107E-03 

.  .  -PHn-4. 

.1.-7.  J 

c 

-2^11_976EL-03__ 

•V 

B-4 

_ly 

.*  "  *  *«•  *  *  *  *  •  'v  '  ■  *  •  ‘  i'*  *  •  .  ,  '  ,  .  '•  ,  *  ,V  •  ,%  J 

r; 

v\1’ 1 T1 .■  TV  *.v J .’7 r*rr\ 

k'lv 

PHI (-4 

‘.-6* 

> 

=  -2. 50772E-03 

PHI (-4 

,-5 

> 

*  -2. 93574E-03 

PHI <-4 

f4 

> 

=  -3.3327E-03 

38 

\-\  PH  I  <  -4 

,-3 

> 

=  -3.52208E-03 

PHI (-4 

,-2 

> 

=  -3.18245E-03 

PHI (-4 

.-1 

> 

=  -1.97734E-03 

PH  I  <  -4 

,  0 

) 

=  -3. 19193E-09 

K-.-; 

f"''.  PHI (-4 

,  1 

> 

=  1 . 97734E-03 

V“-/  PHI (-4 

,  2 

> 

■  3.18245E-03 

AX 

PH I < -4 

.  3 

> 

=  3.52208E-03 

jR 

PHI <-4 

,  4 

> 

=  3. 3327E-03 

PHI (-4 

,  5 

) 

=  2. 93575E-03 

PH  I  <  -4 

i  6 

> 

=  2. 50772E-03 

r; 

PHI <— 4 

.  7 

> 

=  2. 1 1976E-03 

PHI (-4 

,  8 

) 

=  1.79107E-03 

\V- 
o  -. 
•-V 

PHI  <  -4 

,  9 

> 

=  1.5201  IE-03 

PHI <-4 

*  10 

)  =  1 . 29875E-03 

Sv 

PHI < -3 

,-10 

)  =  -1 . 42324E-03 

%xm 

PHI <-3 

,-9 

) 

=  -1 . 69637E-03 

PHI (-3 

> 

=  -2. 04672E-03 

A* 

PHI (-3 

,  “7 

> 

=  -2.49971E-03 

V 

PHI (-3 

» "6 

> 

=  -3.08515E-03 

Vv 

PHI  <  -3 

,  “5 

> 

=  -3. 82592E-03 

PHI <-3 

» -4 

> 

*  -4.6961E-03 

PHI <-3 

,  -3 

> 

=  -5 . 49808E-03 

PHI <-3 

i-2 

) 

=  -5.60857E-03 

PH I < -3 

,-l 

> 

=  -3. 89514E-03 

PHI (-3 

,  0 

> 

=  -6. 59506E-09 

R; 

PHI (-3 

,  1 

> 

*  3 . 895 1 4E-03 

t:v: 

PHI (-3 

»  2 

> 

*  5.60858E-03 

!82 

PHI <-3 

,  3 

> 

=  5.49808E-03 

K* 

PHI (-3 

.  4 

) 

*  4.696UE-03 

^  PHI <-3 

,  5 

> 

*  3. 82592E-03 

sm\  O  phi  <-3 

y\F  PHI  <-3 

*  6 

> 

a  3.08515E-03 

15? 

,  7 

> 

=  2.49971E-03 

k  PHI <-3 

.  8 

> 

a  2. 04672E-03 

K-'-: 

PHI (-3 

,  9 

) 

»  1 . 69637E-03 

tv- 

PH I < -3 

,  10 

)  =  1 . 42324E-03 

!» 

PHI <-2 

,-10 

)  =  -1 . 52488E-03 

•BS 

PHI (-2 

.-9 

> 

*  -1 . 84426E-03 

■ 

PHI <-2 

i-8 

) 

*  -2.26907E-03 

vs 

PHI (-2 

,  -7 

> 

a  -2.8466E-03 

ivC% 

PHI < -2 

t  “6 
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7 
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-1.07131E-03 
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7 
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-8. 13317E-04 
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7 
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8 
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8 
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-8. 95538E-04 
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PHIL-8- 

_9,57962E-_04- 

B-8 

PHI  < 

8  ,  6  > 

= 

9. 69135E-04  ' 

PHI  < 

8,7) 
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=  -6. 73228E-04 
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= 

-7. 47712E-04 
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= 

-7.  68656E-04 
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-7.  41405E-04 
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-6.  75604E-04 

PHI  < 
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= 

2. 1596E-04 
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= 
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PHI  ( 
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PHI  ( 
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7. 6947E-04 

PHI  < 
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= 
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7. 47713E-04 
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PHI  < 
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) 
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) 
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=  6. 16427E-04 
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) 
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PHI  ( 
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) 

=  6.23269E-04 

PHI  < 
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) 

-  6. 05905E-04 

PHI  ( 
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l  ) 

»  5 . 80624E-04 

TABLE  B-l 


Comparison  of  Predicted  Versus  Measured  Surface 

Potential 


Coordinates 

Predicted 

Measured 

Fact 

(x,y) 

mV 

mV 

-10,-10 

0.581- 

90.0 

155. 

-9,-9 

0.714- 

70.0 

98.0 

-8,-8 

0.898- 

90.0 

100. 

-7,-7 

1.165- 

84.0 

72.0 

-6,-6 

• .566- 

114. 

73.0 

-5,-5 

2.211- 

92.0 

42.0 

-4,-4 

3.333- 

144. 

43.0 

-3,-3 

5.498- 

140. 

25.0 

-2,-2 

10.18- 

94.0 

9.0 

-1,-1 

19.06- 

98.0 

5.0 

0,0 

0.000 

00.0 

0.0 

1.1 

19.06 

110. 

6.0 

2,2 

10.18 

112. 

11.0 

3,3 

5.498 

102. 

19.0 

4,4 

3.333 

29.0 

9.0 

5,5 

2.211 

96.0 

43.0 

6,6 

1.566 

110. 

70.0 

7,7 

1.165 

108. 

93.0 
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DOG  AUDITORY  EVOKED  RESPONSE  PLOTS 
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APPENDIX  E 


DOG  AUDITORY  EVOKED  RESPONSES 


On  the  last  day  we  collected  data  from  Ricky,  the 
laboratory  beagle,  we  were  able  to  setup  an  auditory  evoked 
response  (AER)  experiment.  We  used  two  signal  generators  and  a 
pseudo-random  noise  generator  to  generate  a  short  tone  at 
random  interstimulus  intervals.  This  tone  was  coupled  to  the 
dog  with  air  filled  earphones.  The  rest  of  the  setup  was  the 
same  as  for  the  visual  evoked  response  experiment,  that  is, 
the  SQUID  gradiometer  was  positioned  at  different  locations 
over  the  beagles  head  too  get  the  best  response.  The  output  of 
the  SQUID  gradiometer  was  once  again  sent  to  the  Nicolet  660B 
FFT  Analyzer,  so  that,  50  trails  could  be  averaged.  That 
average  was  then  sent  out  of  the  Nicolet  660B  FFT  Analyzer  to 
an  Intel  80/20  computer  so  that  it  could  be  stored  on  a  8  inch 


disk . 


Figure  E-l  shows  an  auditory  evoked  response  taken  from 


a  human  by  electroencephalographic  means,  as  compared  to  an 
auditory  evoked  response  taken  from  a  human  by 
magnetoencephalographic  means,  which  is  shown  in  Figure  E-2. 
In  my  literature  search,  I  was  unable  to  find  an  auditory 
evoked  response  for  a  dog  which  had  been  taken  by 
electroencephalographic  means,  but  was  able  to  find  an 
auditory  evoked  response  for  a  cat  which  had  been  taken  by 
electroencephalographic  means.  This  auditory  evoked  response 
of  the  cat  is  shown  in  Figure  E-3.  The  cat's  brain  is 


physically  smaller  than  the  brain  of  a  beagle,  so  I  would 
expect  the  latencies  measured  from  the  cat  to  be  shorter  than 
the  latencies  we  would  measure  from  the  beagle.  The  plots 
MEG. 066  thru  MEG. 071,  which  are  at  the  end  of  this  appendix, 
shows  the  auditory  evoked  responses  taken  from  the  beagle  by 
magnetoencephalographic  means. 


LATE  AUDITORY  EVOKED"  RESPONSE! 

PROCESSING  IN  PRIMARY  AUDITORY  PROJECTION  AREAS 
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FIGURE  E-l 

HUMAN  ELECTROENCEPHALOHRAPHIC  AUDITORY  EVOKED  RESPONSE 


FIGURE  E-2 


HUMAN  MAGNETOENCEPHALOGRAPHIC  AUDITORY  EVOKED  RESPONSE 


FIGURE  E-3 

CAT'S  ELECTROENCEPHALOGRAPHIC  AUDITORY  EVOKED  RESPONSE 

In  the  case  of  the  human  auditory  evoked  response  data, 
we  see  that  both  the  electroencephalographic  auditory  evoked 
response  and  the  magnetoencephalographic  auditory  evoked 
response  have  peaks  at  latencies  of  approximately  50 
milliseconds  and  180  milliseconds. 

For  the  cat  versus  the  beagle  case,  we  first  see  that 
the  cat  has  a  negative  going  waveform  during  the  first  20 
milliseconds  of  its  response.  As  expected,  the  beagles 
waveform  is  going  negative  out  to  40  milliseconds. 
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